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airborne  computers  used  in  Air  Force  avionic  weapon  systems.  The 
routines  implemented  include  several  algebraic  functions  that  are  intended 
to  serve  as  a  benchmark  for  future  contractor  development.  Appendix  A 
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Abstract 

This  project  produced  a  run-time  math  library  for  the  MIL-STD-1750A 
embedded  computer  architectures.  The  math  library  consists  of  the 
algebraic  functions.  In  addition,  the  steps  required  for  a  performance 
analysis  of  the  math  library  have  been  outlined. 

Several  approximation  methods  were  investigated.  The  Chebyshev 
Economization  of  Maclaurin  series  polynomials,  and  rational 
approximations  derived  from  the  second  algorithm  of  Remes  were 
determined  to  be  the  best  methods  available.  Each  function’s 
implementaion  was  designed  to  take  advantage  of  features  of 
M1L-STD-1750A  architectures.  The  recommended  test  procedures  provide 
measures  of  the  average  and  worst  case  generated  errors  within  each 
approximation. 


The  Air  Force  is  interested  in  reducing  the  life-cycle  costs  of  its 
avionics  weapon  systems.  Standardization  of  high  order  languages  and 
an  Instruction  Set  Architecture  (ISA)  are  two  of  the  many  ways  the  Air 
Force  can  reduce  these  costs.  In  the  past,  a  major  cost  contributor  was 
the  proliferation  of  unique  avionics  systems  and  subsystems.  Costs 
increased  with  respect  to  purchasing  and  inventorying  small-lot  spares 
at  many  bases,  training  technicians  to  maintain  complex  and/or  unique 
flight  and  test  equipment,  developing  and  maintaining  software 
development  facilities,  training  programmers  to  write  application 
programs  in  seldom  used  high  order  languages,  and  training  programmers 
to  maintain  software  (especially  operating  systems)  in  seldom  used 
machine  languages.  ( 1  -  8.1) 

M1L-STD-I750A  defines  a  standard  16-bit  instruction  set 
architecture  intended  primarily  for  avionics  weapon  systems.  The  major 
cost  advantage  of  this  standard  ISA  comes  in  the  form  of  common 
support  software  tools.  An  extensive  set  of  support  software  tools  has 
already  been  developed  and  includes  a  1750A  assembler /crossassembler, 
a  J73  compiler  with  I750A  ISA  code  generator,  a  linker  program,  a 
loader  program,  and  a  1750A  acceptance  test  program.  (I;  8.4)  Other 
cost  benefits  are  realized  through  the  independent  development  of 
software  and  hardware,  (2:  1)  and  common  maintenance  and  test 
equipment.  (3:  168) 


Standardization  of  languages  also  has  an  impact  on  cost  reduction. 
"In  1978,  the  Department  of  Defense  had  in  its  inventory,  software 
written  in  about  150  different  programming  languages.  This  linguistic 
proliferation  increased  maintenance  problems  due  to  programmer  training 
requirements  and  lack  of  support  tools  for  many  of  the  languages." 
(1-  6.1)  The  D.O.D.  and  Air  Force  recognized  this  as  a  problem,  and  they 
took  steps  to  correct  it.  The  D.O.D.  instruction  5000.31,  "Interim  List  of 
D.O.D.  Approved  High  Order  Programming  Languages,"  states  that  only 
approved  languages  may  be  used  for  new  defense  system  software. 
JOVIAL  is  one  of  the  few  languages  approved  by  this  instruction. 

As  previously  mentioned,  the  development  of  a  standard  ISA  such  as 
MIL-5TD-1750A  helps  reduce  total  life-cycle  costs  for  Air  Force 
avionics  weapon  systems.  This  reduction  is  partially  due  to  the  use  of 
common  support  software  tools,  many  of  which  have  already  been 
developed.  As  was  previously  mentioned,  one  of  the  support  software 
tools  that  has  been  developed  is  a  JOVIAL  compiler  that  generates 
1750A  ISA  code.  However,  a  math  library  containing  the  algebraic  and 
trigonometric  functions  required  by  this  language  has  not  been  developed. 
The  sponsor  for  this  thesis  is  the  Aeronautics  System  Division,  Language 
Control  Branch.  They  are  the  D.O.D.  JOVIAL  and  ADA  compiler  validation 
site,  and  are  responsible  for  the  development  of  such  libraries. 
Completion  of  this  thesis,  with  its  development  of  a  math  library  for 
software  support  of  1750A  systems,  can  help  the  Air  Force  reduce 
avionics  weapon  systems  costs. 


Problem 


Prior  to  the  completion  of  this  thesis  effort,  there  were  no  math 
libraries  written  to  take  advantage  of  the  1750A  instruction  set.  In 
keeping  with  the  intent  of  recent  standardization  policies  of  both  the 
D.O.D  and  Air  Force,  the  library  created  by  this  thesis  is  written  in  the 
D.O.D  approved  language  JOVIAL.  Actual  coding  of  the  library  was  only  a 
small  part  of  this  thesis.  Most  of  the  detail  has  gone  into  verification, 
validation,  and  evaluation  of  the  product.  As  such,  the  focus  of  this 
report  is  divided  into  two  primary  categories:  software  development  and 
software  testing. 

Math  libraries  are  important  because  they  provide  the  programmer 
several  tools  that  serve  as  building  blocks  for  applications.  Math 
libraries  prevent  programmers  from  having  to  recreate  each  function 
whenever  one  is  required  for  use.  Libraries  also  provides  a  means  for 
using  functions  that  take  full  advantage  of  the  computer  architecture  for 
which  they  were  written. 

The  design  of  a  procedure  for  computing  the  value  of  functions  is 
not  mathematically  complete  unto  itself.  An  understanding  of  a 
computer  architecture's  operation  is  necessary  to  insure  that  the 
computation  of  any  given  function  is  as  efficient  as  possible,  while  also 
providing  the  highest  degree  of  accuracy.  Such  architectural 
considerations  include  word  size,  number  of  bits  in  both  the  exponent 
and  coefficient  fields  of  a  floating  point  number,  the  way  mathematical 
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operations  are  performed  by  the  architecture,  memory  size  of  the 
architecture,  and  execution  time.  Other  considerations  include  overflow, 
underflow,  and  precision.  These  considerations  for  the  mathematic 
functions  define  some  of  the  problems  addressed  by  this  thesis  effort. 

Scgpg 


This  effort  was  limited  to  the  design,  code,  and  evaluation  of 
algebraic  functions.  The  functions  were  included  in  a  math  library 
targeted  for  MIL-STD-I750A  computer  architectures,  and  are  the  ones 
typically  found  in  most  FORTRAN  libraries.  Specifically,  these  functions 
include  square  root  (sqrt),  exponential  (exp),  natural  logarithm  (alog), 
and  common  logarithm  (alog  10). 

All  functions  have  been  written  to  accept  and  return  only  extended 
precision  floating-point  values.  The  specific  floating-point  functions 
are  invoked  by  using  the  name  given  above. 

Performance  summaries  for  each  of  the  functions,  and  algorithms 
are  provided  in  Appendix  E.  They  may  be  used  to  determine  the 
polynomial  coefficients  for  computing  any  of  functions  addressed  by  this 
paper.  These  algorithms  produce  coefficients  that  are  valid  for  any 
nonvector  architecture. 
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During  a  design  review  held  in  nay  of  1985.  it  was  made  clear  that 
certain  events  could  cause  overflow  and  underflow  errors,  and  division 
by  zero.  Since  the  functions  are  to  be  used  within  an  embedded  avionics 
weapon  system,  it  is  necessary  that  such  conditions  are  detected  and 
handled  gracefully.  The  consensus  of  opinion  from  all  participants  of  the 
design  review  was  that  the  functions  should  not  be  aborted,  and  that 
default  values  should  be  returned.  The  error  conditions  and  values 
returned  are  discussed  in  the  individual  design  sections  of  this  thesis. 
This  constitutes  an  important  assumption  on  how  to  handle  such  error 
conditions,  and  needs  further  investigation  before  implementation  on  a 
real-time  system. 

Another  factor  discussed  during  the  design  review  was  the 
distinction  between  fixed-point  and  floating-point  functions. 
Floating-point  functions  have  greater  precision  than  fixed-point 
algorithms,  but  take  longer  to  execute.  Although  the  fixed-point 
functions  are  faster,  the  algebraic  math  routines  and  the  JOVIAL 
computer  language  do  not  lend  themselves  to  this  method  of  calculation. 
Therefore,  as  stated  earlier,  only  the  floating-point  algebraic  math 
library  functions  have  been  implemented. 

Qgnecal  AgpcoasH 

The  approach  used  during  this  thesis  effort,  is  termed  the 
"logical ized"  model  of  a  software  system  development  cycle.  This 


approach  was  considered  a  better  alternative  than  the  more  commonly 
used  "waterfall"  method  of  software  system  development.  The 
“waterfall"  method  is  composed  of  a  neat,  concise  and  logical  ordering  of 
a  series  of  steps,  each  of  which  must  be  accomplished  in  order  to  obtain 
a  final  software  product.  These  steps  are  performed  in  order  and  include 
systems  analysis,  requirements  definition,  preliminary  design,  detailed 
design,  coding,  testing,  and  implementation. 

The  “logicalized"  model  is  similar  to  the  “waterfall"  model  just 
described,  but  it  is  more  concerned  with  the  problem  definition  part  of 
the  cycle  (see  Figure  1).  This  approach  is  more  useful  in  eliminating 
errors  that  typically  occur  during  the  requirements  definition  and  design 
phases  of  the  “waterfall"  method.  Errors  generated  during  these  phases 
typically  occur  because  designers  have  a  tendency  to  shift  between 
abstract  high-level  design  issues  and  physical  implementation 
considerations.  Thayer  (5;  335-41)  and  Boehm  et  al.  (6:  125-33) 
made  it  clear  that  these  problems  exist,  and  that  design  errors  not  only 
outnumber  other  errors,  but  that  they  are  also  more  persistent.  For  this 
reason,  more  attention  was  given  to  the  top-down  decomposition  and 
abstract  (logical)  modeling  of  this  particular  software  system.  Such  a 
structured  approach  recommends  a  dichotomy  between  the  logical  design 
issues,  and  implementation  issues. 
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Figure  1  Chart  of;  (a)  Waterfall;  and  (b)  Logicalized 
Model  of  a  Software  Development  Cycle 
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Table  1  Information  Flow  of  the  Logical ized  Software 
Development  Cycle 


The  information  flow  of  a  logical  ized  model  is  summarized  in 
Table  1.  and  is  “analogous  to  an  artist's  conception  of  a  building,  i.e. 
there  is  enough  information  to  allow  the  customer  and  designer  to 
communicate  and  to  establish  the  buildings  pluses  and  minuses,  but  not 
enough  detail  to  begin  construction.  A  series  of  reviews,  refinements, 
and  the  imposition  of  local  building  ordinances,  for  example,  are 
necessary  before  construction  can  start."  (7=14) 

Therefore,  the  approach  taken  for  this  project  was  similar  to  that 
just  described.  The  ASD/Language  Control  Branch  established  the 
requirements  for  a  MIL-STD-1750A  run-time  math  library  written  in  the 
D.O.D  approved  high-order  language  JOVIAL.  During  a  design  review  and 
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several  other  meetings,  certain  design  considerations  were  refined. 
Then  a  "logical”  model  was  established  as  a  baseline.  This  was 
accomplished  by  using  the  refined  problem  statement,  and  researching 
the  different  methods  for  approximating  the  different  algebraic 
functions. 

The  baseline  model  served  as  a  reference  from  which  all  decisions 
regarding  actual  implementation  could  be  made.  Before  proceding  to  the 
next  phase  of  development,  two  such  decisions  had  to  be  made.  These 
decisions  were  to  determine  which  testing  methods  and  which 
performance  evaluation  techniques  would  be  used  after  coding  was 
complete.  These  decisions  determ ined  what  sort  of  tests  would  catch 
all  possible  errors,  and  determined  what  techniques  could  be  used  to 
establish  a  confidence  level  for  the  final  product. 

Up  until  this  point,  the  abstract  model  has  been  devoid  of  any 
implementation  considerations.  However,  after  it  was  clear  that  the 
abstract  design  was  complete  and  consistent  with  the  requirements,  it 
became  necessary  to  consider  changes  to  fit  the  problem  into  the 
n IL-STD- 1 750A  environment.  Before  any  changes  could  be  made,  it  was 
necessary  to  complete  the  following  steps:  study  the  architecture  and 
ISA  defined  by  niL-STD-1750A;  determine  what  resources  were 
available,  such  as  software  support  tools  and  hardware;  and  then  to  learn 
how  to  use  the  available  resources.  From  there,  it  was  possible  to 
develop  an  abstract  model  of  an  implementable  solution.  This  model 
took  advantage  of  those  environmental  factors  that  affected  the  speed 
and  accuracy  of  computation  for  each  function  approximation. 
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The  major  subset  of  the  logical ized  software  engineering 
methodology  just  described  is  called  structured  programming. 
Structured  programming  can  be  understood  as  the  decomposition  of  a 
problem  in  order  to  establish  a  manageable  problem  structure.  The 
highest  conceptual  level  represents  a  general  description  of  the  problem, 
and  each  level  of  decomposition  provides  more  detail  into  the  problem. 
This  decomposition  is  carried  out  until  the  problem  is  almost  in  coded 
form,  and  is  often  called  a  stepwise  refinement  of  the  problem.  All 
implementation  considerations  are  left  until  the  lowest  levels  of 
refinement. 

The  goals  of  structured  programming  must  be  to  minimize:  the 
number  of  errors  that  occur  during  the  development  process;  the  effort 
required  to  correct  errors  in  sections  of  code  found  to  be  deficient;  the 
effort  required  to  upgrade  sections  when  more  reliable,  functional,  or 
efficient  techniques  are  discovered;  and  the  life-cycle  costs  of  the 
software.  (8=  32)  It  must  also  reduce  the  complexity  of  the  problem. 

Structured  flowcharting  is  a  technique  used  to  support  these 
structured  programming  concepts  and  goals,  and  is  “designed  to  reduce 
labels  and  unstructured  branching,  encourage  a  single  entry/single  exit 
approach,  aid  in  the  use  of  top-down  design  techniques,  and  enhance 
modularization.  The  approach  encourages  the  designer  to  conceive  of  the 
system  in  high-level  constructs  and  not  in  terms  of  individual  detailed 
statements."  (7:  116)  The  structured  flowcharting  technique  was  used 
throughout  the  development  of  this  project,  not  only  because  of  the 
reasons  just  mentioned,  but  also  for  its  simplicity  and  understandability 
from  a  reviewer  standpoint. 
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This  thesis  addresses  the  design  and  performance  evaluation  of  a 
run-time  math  library  that  is  targeted  to  I1IL-STD-1750A  architectures. 
The  requirements  definition  for  this  problem  has  already  been  discussed 
(Chapter  l  -  Problem/Scope).  The  next  topic  discussed  is  the  theoretical 
development  of  this  thesis  effort  (Chapter  2).  This  is  followed  by  a 
discussion  of  the  detailed  design  considerations  that  were  made  during 
implementation  of  the  library  functions  (Chapter  3).  The  last  aspects 
covered  in  this  report  are  the  test  and  performance  evaluation  methods 
considered  (Chapter  4)  and  conclusions  (Chapter  5). 

Appendices  include  algorithms  useful  for  determining  the 
pseudo-code  operations  used  in  the  structured  flowcharts  (Appendix  A), 
source  listings  for  the  implemented  functions  (Appendix  B),  support 
software  developed  in  conjunction  with  this  thesis  (Appendix  C),  the  VAX 
VMS  command  files  required  to  compile,  link,  and  run  the  developed 
product  (Appendix  D),  and  the  coefficients  for  each  of  the  functions 
(Appendix  E). 


The  purpose  of  this  thesis  was  to  create  and  analyze  algebraic 
functions  developed  for  1750A  architectures.  This  chapter  is  concerned 
with  the  design  theory  of  the  algorithms  used  to  approximate  those 
functions.  Within  the  given  constraints,  the  emphasis  for  each  of  the 
designs  is  to  compute  results  as  quickly  and  as  accurately  as  possible. 

One  way  of  computing  a  value  quickly  is  to  select  an  approximation 
that  converges  rapidly  towards  the  value  of  the  true  function,  f(x).  There 
were  several  methods  of  approximation  that  were  considered;  however,  the 
polynomial  and  relational  approximations  described  by  Cody  and  Waite 
(4--  17-84)  were  found  to  be  the  best.  The  coefficients  given  by  Cody  and 
Waite  were  derived  by  using  Chebyshev  Economization  of  the  Taylor  series 
for  each  function  for  the  approximation  itself,  or  as  a  starting  point  for 
computing  a  rational  approximation  via  the  second  algorithm  of  Remes. 
An  excellent  reference  for  Chebyshev  Economization  is  Conte  and  de  Boor 
(9;  265-273),  and  an  excellent  reference  for  the  second  algorithm  of 
Remes  is  Ralston.  ( 1 0*  301-306) 

Another  means  of  reducing  the  amount  of  processing  time  required  to 
compute  a  result  is  to  take  advantage  of  certain  aspects  of  the  computer’s 
architecture,  as  well  as  the  different  execution  times  for  different 
instructions  within  the  ISA.  For  example,  incrementing  the  exponent  field 
of  a  floating-point  value  is  not  only  faster,  but  more  accurate  than  the 


equivalent  operation  of  multiplying  by  two,  or  examining  the  sign  bit  of  a 
variable  is  faster  than  comparing  the  entire  value  to  zero.  These 
techniques  have  been  used,  and  are  referenced  in  the  design  descriptions 
as  pseudo-operations.  These  operations  are  equivalent  to  those  described 
by  Cody  and  Waite  (-4-'  9),  and  are  listed  in  Appendix  A. 

The  accuracy  of  an  approximation  may  be  dependent  upon  the  domain 
over  which  the  function  is  approximated.  For  example,  if  the  domain  of  an 
approximation  is  halved,  the  error  may  be  reduced  by  a  factor  of  about 
2-<ft+i>  for  au  polynomials  of  degree  n.  (1 1=  59)  This  can  be  shown  to  be 
true  for  most  functions,  but  not  all  of  them.  Domain  reduction  has  no 
effect  on  accuracy  in  approximations  of  certain  functions;  however,  it 
still  serves  as  an  excellent  guide  when  designing  an  application.  This  is 
due  to  the  way  computer  architectures  perform  operations  and  store 
mathematical  values  for  floating-point  numbers.  The  most  significant 
bits  of  a  number  are  always  maintained,  and  since  only  a  finite  number  of 
bits  are  available  to  represent  the  value,  it  is  possible  that  bits  from  a 
fractional  representation  may  be  lost  during  operations  on  large  numbers. 

Approximation  Techniques 


The  MIL-STD- 1 750A  ISA  doesn't  call  for  the  implementation  of  the 
elementary  functions  as  standard  instruction  operators,  so  it  is  necessary 
to  design  software  routines  of  optimum  efficiency  to  replace  them.  The 
word  "optimum”  could  be  given  a  variety  of  precise  definitions,  but 
presumably  it  refers  to  an  average  execution  time  and  storage  space. 


Unfortunately,  there  is  no  known  way  to  derive  or  prove  such  an  "optimal" 
design.  For  these  reasons,  the  search  for  the  appropriate  approximation 
technique  was  limited  to  polynomial  and  rational  approximations. 

Some  of  the  most  popular  methods  of  approximation  used  are  called 
CheDyshev  approximations.  Chebyshev  approximations  are  often  referred 
to  as  "minimax"  approximations  because  they  are  used  to  minimize  the 
maximum  "error"  between  the  true  function  f(%),  and  the  approximation 
of  f(x).  However,  these  methods  of  approximation  are  not  without  their 
problems,  and  there  is  a  price,  even  though  it  is  small,  to  be  paid  for 
using  them.  For  example,  the  sum-of-squares  of  the  errors  in  a  Chebyshev 
approximation  will  be  higher  than  if  a  least-squares  method  of 
approximation  is  used.  However,  since  Chebyshev  approximations  assure 
that  an  error  is  never  greater  than  a  given  amount,  they  were  selected  by 
this  study. 

Polynomials.  The  first  class  of  approximations  discussed  are 
polynomials,  and  are  the  simplest  of  all  the  classes  of  approximations 
considered.  The  most  important  subclass  of  the  polynomials  is  the  class 
(Chebyshev),  and  are  polynomials  not  exceeding  degree  n.  The 
Chebyshev  polynomials  are  especially  important,  and  gave  rise  to  the 
general  concept  of  Chebyshev  "approximations"  discussed  in  the  preceding 
paragraph. 

The  motivation  for  using  Chebyshev  polynomials  over  all  other 
polynomials  is  their  property  of  least  maximum  error,  and  their  error 
behavior  over  the  entire  interval  of  the  approximated  function.  Through 
the  use  of  Theorem  I.  the  Alternation  Theorem  given  below,  Chebyshev 
was  able  to  prove  for  all  the  polynomials  of  degree  n  with  a  leading 


coefficient  of  1,  that  the  Chebyshev  polynomial  divided  by  2  has  the 
least  maximum  error  in  the  interval  [-1,1].  In  other  words,  no  other 
polynomial  of  the  type  mentioned  will  have  a  smaller  error  than 
Tft(x)/2ar1.  In  order  for  a  polynomial  Pft(%)  to  be  considered  a  Chebyshev 
approximation  of  the  function  f(x),  the  theorem  requires  that  the 
maximum  discrepancy  between  j(x)  and  Pft(x)  occur  with  alternating 
signs  at  n+2  points  over  the  interval  [-1,1]. 

Alternation  Theorem:  The  polynomial  Pn  of  degree  in  that  (I) 

best  approximates  f  is  characterized  by  the  existence  of  at 
least  n*2  "points  of  alternation" 

The  other  motivation  for  the  use  of  Chebyshev  polynomials  is  that 
its  generated  errors  are  more  well  behaved  than  the  errors  generated  by 
other  polynomials.  For  example,  approximations,  based  on  the  Maclaurin 
series  whose  interval  includes  zero,  have  errors  that  are  very  nonuniform 
--  small  near  the  middle,  but  very  large  at  the  end  points.  It  is  more 
desirable  to  use  an  approximation  whose  behavior  is  more  uniform  instead 
of  powers  of  x.  Since,  as  stated  in  Theorem  1,  the  Chebyshev 
polynomials  spread  the  error  over  the  entire  interval,  they  provide  this 
more  desirable  behavior. 

Definition  qL  Chebyshev  Polynomials.  The  Chebyshev  polynomials 
form  an  orthogonal  set,  and  are  defined  by  the  following  equation. 

T.  (x)  =  cos  (n8)  9  =  arccos  (x)  ( 1 ) 

n  =  1,  2,  .  .  . 
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From  elementary  trigonometry,  cos(n8)  is  a  polynomial  of  degree 
n  in  cos(e),  and  ccs(arccos(x))  =  %;  therefore,  it  follows  that  the 
Chebyshev  polynomials  defined  by  xft(x)  =  cos(n  arccos(x))  are 
polynomials  of  degree  n. 

By  substituting  arccos(x)  for  0  and  xft(x)  for  cos(n  arccos(x)) 

in  the  identity  function  shown  in  equation  (2),  the  recurrence  relation 
defined  in  (3)  is  formed. 

cos  ((n  ♦  1)0)  ♦  cos  C(n  -  1)0)  =  2  cos  (0)  cos  (n8)  (2) 

=  2x  xm(x)  -  xa_,(x)  (3) 

Let  x0  =  1  and  x,  =  x.  then  from  the  recurrence  relation 

defined  in  (3).  successive  polynomials  of  greater  degree  can  be  generated 
as  in  column  A  of  Table  2. 

By  using  the  results  in  column  A  of  Table  2,  the  powers  of  the 
Chebyshev  polynomials  can  be  found.  That  is,  it  is  possible  to  express 
the  powers  of  x  in  terms  of  xft.  An  example  of  the  powers  of  x  are 

shown  in  Table  2  column  ft  Appendix  A  contains  an  algorithm  that 
generates  both  the  Chebyshev  polynomials,  and  their  powers. 


1  .r'"  tT9  yi 


A 

£ 

To  =  1 

T,  =  x 

x°  = 

X1  = 

1  =  T0 

x  =  T, 

T2  =  2% 

Ti  -  T0  -  2x*  -  l 

X"  = 

Vx  (2x‘t  -  1  +  1  )  =  V2  ( x  2  +  xn) 

r3  r  2x 

r2  -  T,  r  4x-  -  3x 

X3  = 

(4x*  -  3x  +  3x  )  =  %(x7  +  3t.) 

Table  2  (A)  Chebyshev  Polynomials;  (B)  Powers  of  Chebyshev  Polynomials 


Chebyshev  Economization.  As  already  mentioned,  the  flaclaurin 
series  can  be  used  to  approximate  many  functions.  In  addition  to  the 
disadvantages  that  have  already  been  mentioned  for  using  this  series  as  an 
approximation,  the  Maclaurin  series  also  converges  very  slowly.  That  is, 
it  takes  several  multiplications  and  additions  to  obtain  a  desired 
accuracy.  One  way  of  obtaining  a  lower  degree  polynomial,  and  still 
maintain  the  desired  accuracy,  is  to  use  a  technique  that  is  called 
■'telescoping"  or  “Chebyshev  Economization".  In  other  words,  the 
polynomial  can  be  expressed  in  a  manner  similar  to  that  shown  in  (4). 

Pn(x)  =  d0x0(x)  +  .  .  .  +  daxa  (4) 


To  compute  the  economized  polynomial  approximation  to  the  function 
f(x)  of  absolute  accuracy  e  on  the  interval  [-1,  I],  use  the  following 
procedure  as  outlined  by  Conte  et.al.  (9:  271-272) 


17 


Step  1.  Get  a  power  series  expansion  for  j(x)  valid  on  [-1,1]; 
typically,  calculate  the  Maclaurin  or  Taylor  series  expansion  for  f(x) 
around  x  =  0. 

Step  2.  Truncate  the  power  series  to  obtain  a  polynomial  as  in  (5). 
which  approximates  f(x)  on  [-1,  !]  within  an  error  where  is 

smaller  than  e,  and  e.  is  defined  as  in  (6).  The  result  of  e.  is  the 
maximum  absolute  value,  within  the  interval  [-1,1],  of  the  product  of 
the  first  truncated  coefficient,  x  to  the  power  of  n  +  1,  and  the  n  +  1 
derivative  of  the  function  f(x). 

Pa(x)  -  a0  +  atx  +  .  .  .  +  aftxtt  (5) 

=  =  f  (a*1J(x)  (6) 

Step  3-  By  making  use  of  a  table  similar  to  that  shown  in  Table  2 
column  a,  expand  the  polynomial  Pfl(x)  into  a  Chebyshev  series  as 

defined  in  (4).  In  other  words,  substitute  the  far  right-hand-side  of  the 
equations  in  Table  2  column  B,  with  the  appropriate  powers  of  x 
contained  in  the  polynomial  formed  by  Step  2  of  this  algorithm.  The 
result  is  similar  to  that  shown  in  (7),  but  of  a  greater  degree. 


.1  M  >.l 


V K'  ■. 


Step  4.  Retain  the  first  k  +  1  terms  in  this  series,  i.e.  find 
equation  (7),  choosing  k  as  the  smallest  possible  integer  such  that 
equation  (8)  holds  true. 

P  k(x)  =  d0x0(%) +  .  .  .  +  d^Tjj  (7) 

e.  *  '  •  ■  ■  *  d,  5  E  (8> 

Step  5.  Convert  the  result  of  Step  4  into  a.  power  series  polynomial 
similar  to  (5),  by  making  use  of  a  table  similar  to  that  in  Table  2  column 
A.  In  other  words,  substitute  the  right-hand-side  values  of  Table  2 
column  A,  into  the  equation  formed  by  Step  4.  Simplify  the  result. 

Rational  Approximations.  In  most  instances,  rational  approximations 
will  generate  a  least  maximum  error  that  is  as  small  or  smaller  than  a 
Chebyshev  polynomial,  and  will  also  cost  less  in  terms  of  the  number  of 
multiplications  and  additions  required  to  compute  them.  Therefore,  they 
deserved  attention  in  this  study. 

As  stated  earlier,  the  approximation  techniques  considered  by  this 
thesis  are  classified  as  Chebyshev  approximations.  These  methods, 
through  their  exploitation  of  Theorem  1,  provide  approximations  whose 
maximum  error  is  less  than  those  generated  by  other  techniques.  There 
are  several  algorithms  that  generate  rational  approximations  that  can  be 
considered  Chebyshev  approximations;  however,  the  ones  that  generate  the 
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most  uniform  approximations  are  those  generated  by  the  second  algorithm 
of  Remes.  This  algorithm  is  easily  automated,  and  is  described  in  detail 
by  the  following  subsection. 

]M  Second  Algorithm  qL  Remes.  The  method  used  in  this 
description  is  similar  to  that  outlined  by  Ralston  (10:  301-305),  and  is 
summarized  in  Figure  2. 

Let  }(x)  be  a  continuous  function  that  is  to  be  approximated  over  the 
the  interval  [a,  b],  and  let  the  interval  include  the  point  0.0. 
Furthermore,  let  (9)  equal  the  error  of  any  rational  approximation  of  the 
form  shown  in  (10). 


rmjk(x)  =  max  *  ft**  "  V*(x}  I 


(9) 


‘inrk  = 


P  Px'1 
Qk!x) 


m 

Z  Si  x 


i=o 


X  i 
1  b.x 

j=o  > 


( 1 0) 


Step  1  of  the  algorithm  names  the  input  required  for  this  algorithm. 
The  input  value  f(%)  is  the  function  being  approximated.  If  the  algorithm 
is  being  run  on  a  machine  with  higher  precision  than  the  error  for  which 
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the  function  is  being  approximated,  then  the  built-in  functions  of  the 
machine  can  be  used  for  f(%).  If  the  machine  that  the  algorithm  is  to  run 
on  is  of  the  same  precision  for  which  the  approximation  is  to  be  made, 
then  a  reasonable  substitute,  such  as  a  truncated  power  series  that  is  of 
equal  or  greater  precision  than  what  is  being  approximated,  can  be  used. 

The  other  inputs  include1  m,  k,  [a,  b],  and  C0  .  .  . (^ .  The  values 

m,  and  k  represent  the  degree  of  the  polynomials  found  in  the 
numerator  and  denominator,  respectively.  The  interval  [  a,  b  ]  is  the 
interval  for  which  the  approximation  is  valid,  and  should  include  the  point 
zero,  as  it  will  allow  the  coefficient  b0.  of  the  denominator,  to  always 
be  one.  The  values  C0  .  .  .  represent  the  first  N  +  1  coefficients  of 
the  power  series  polynomial  that  is  being  converted  to  a  rational 
approximation.  The  value  N  represents  the  sum  of  the  degree  of  the 
polynomials  used  in  the  numerator  and  the  denominator  (m  +  k). 

The  second  step  of  the  algorithm  is  to  compute  a  series  of  Pade 
approximations  and  their  error  coefficients.  The  Pade  approximations  are 
of  the  form  depicted  in  (11),  with  the  restrictions  that  0  i  i  s  m  and 
0  i  j  -  i  i  k.  For  example,  the  sequence  of  Pade  approximations  computed 
for  an  R22  approximation  would  only  include  R<0>00(x),  R<1>1<0(%), 

R<3>u(x),  and  error  of  the  approximations  is 

equal  to  the  first  power  of  %  truncated  from  the  power  series,  multiplied 
by  the  error  coefficients  shown  in  (12).  The  error  calculations  used  would 
only  include-  d*0*0^,  d*1'0^,  d<u>3,  d*2*1^,  and  d*2-2^  . 


The  coefficients  for  each  of  the  sequence  of  Pade  approximations  are 
computed  using  (13)  and  (14).  Equation  (13)  forms  a  set  of  m  linear 
equations,  which  when  solved,  determines  the  value  of  the  coefficients 
used  in  the  denominator.  Those  values  can  then  be  directly  substituted 
into  the  set  of  equations  formed  by  (14).  and  will  determine  the  value  of 
the  coefficients  for  the  numerator. 

j? 

,?0CN-s-jbi  =  0  s  =  0,  1 , . . . ,  N  -  m  -  1  (13) 

( Cj  -  0  if  j  <  0.,  b0  =  1 ) 

a  =  f  r  b-  r  =  0,1 - m  (14) 

’  (bj=0  if  j  >  k  ) 

The  third  step  of  the  Remes  algorithm  is  to  compute  the  economized 
approximation  <w*>-  To  complete  this  step,  it  is  necessary  to  compute 
the  Chebyshev  polynomial  x^.  This  polynomial  can  be  determined  by 
using  equation  (3)  of  the  previous  subsection.  Once  the  coefficients  of 
Ttbi  are  foun<1-  then  values  y  from  (15)  can  be  directly  substituted 
into  (16),  and  thus  solve  w  %).  The  value  tj  in  (15)  is  the  coefficient 


23 


for  u*  in  t^u).  The  rational  approximation  must  also  be  normalized, 
that  is,  the  numerator  and  denominator  must  be  divided  by  b0,  such  that 
b0  will  remain  equal  to  1. 
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The  final  step  of  the  Remes  algorithm  is  an  iterative  one.  Now  that 
the  initial  approximation  to  the  function  has  been  found,  it  becomes 
necessary  to  find  the  N  ♦  2  points  of  alternation.  This  can  be  done 
through  interpolation,  or  by  dividing  the  interval  into  several  small  pieces 
and  solving  for  each  point  on  a  division.  This  method  works,  and  all  that 
is  necessary  is  a  little  bookkeeping  to  maintain  a  list  of  the  N  ♦  2 
points  of  alternation.  This  step  consists  of  the  following  three 
procedures. 


Procedure  I.  Solve  the  system  of  N  +  2  equations  for  the  N  ♦  2 

unknowns  a0<0> . . bk<°>  ,  and  as  shown  in 

expression  (17).  Note  that  E*0>  is  the  magnitude  of  error  in  the 
approximation  at  each  of  the  points  .  and  for  the  first  iteration 

can  he  assumed  to  be  0. 


fix?0] - =  H)4E  (17) 

A  Vxflj 

Procedure  2.  Find  h0(x)  as  shown  in  ( 18).  The  function  h0(x)  then 

I  •  has  a  magnitude  of  |E**|  with  alternating  signs  at  ,  i=0 . N  *  ). 

In  the  neighborhood  of  each  x^®,  there  is  a  point  x^1*  at  which  h0(x) 
has  an  extremum  of  the  same  sign  as  that  of  f(x)  -  K^^Cx) 
at  x^0>  .  Replace  each  x^0*  by  the  corresponding  x^1*  .  If  x  .the 
point  at  which  h0(x)  has  its  maximum  magnitude,  is  one  of  the  points 
XjW  .  do  not  perform  procedure  3.  If  not,  replace  one  of  the  points  xJV 
by  x  in  such  a  way  that  h0(x)  still  alternates  in  sign  on  the  points 


h0(x)  =  t(x) 


1=0 


x’ 


k 

i  =  o 


X> 


1  ] 


(18) 
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Procedure  3.  Repeat  procedures  1  and  2  using  the  points 

XjW . <n  This  process  generates  a  sequence  of  rational 

approximations  which  will  converge  to  an  optimum  if  the  initial  extrema 
were  sufficiently  close. 


111.  Development  and  Design  of  the  F unci. i ops 

General  Discussion 

This  chapter  deals  with  the  detailed  design  of  each  of  the  specific 
functions.  Each  design  has  an  associated  structured  flowchart,  and  each 
Pox  within  the  flowchart  has  been  numbered  for  ease  of  reference. 
Pseudo-operations  are  used  throughout  each  of  the  flowcharts,  and 
includes  those  defined  by  Cody  and  Waite  (4;  9-10).  Furthermore,  a  few 
additional  pseudo-operations  have  been  introduced,  (see  Appendix  A) 
Although  the  approximation  methods  used  are  those  suggested  by 
Cody  and  Waite,  the  actual  design  implementations  are  significantly 
different.  The  designs  proposed  by  Cody  and  Waite  are  guidelines  for  a 
broad  class  of  computer,  and  weren't  specifically  targeted  towards  a 
1750A  architecture.  Therefore,  the  designs  have  been  tailored  somewhat. 

The  coefficients  for  each  of  the  functions  were  either  taken  from 
Cody  and  Waite,  or  are  modifications  of  those  provided  by  Cody  and 
Waite.  These  modifications  are  discussed  in  their  appropriate  subsection. 

Square  Root  Implementation 


The  square  root  of  every  non-negative  floating  point  number  "X"  can 
be  computed.  Computation  is  composed  of  three  steps:  the  reduction  of 
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the  given  argument  X  into  the  parameters  f  and  “e"  using  base  2, 


X  =  f  *  2e, 


1/2  i  f<l 


sqrt(X)  -  sqrt(f)  *  2e/^, 


if  e  is  even 


sqrt(X)  -  [  sqrt(f)  /  sqrt(2)  1  *  if  e  is  odd 


the  computation  of  sqrt  (f),  and  the  reconstruction  of  sqrt  (X)  from  the 
results. 

The  variable  “X"  is  the  argument  passed  to  the  square  root  function. 
Since  JOVIAL  treats  formal  arguments  as  read  only,  upon  program  entry, 
the  value  of  “X"  is  assigned  to  T*.  (step  1  of  Figure  3)  The  variable  "  F" 
is  then  used  throughout  the  remainder  of  the  procedure. 

The  next  step  is  to  check  if  ”F“  is  either  zero  or  one.  (step  2  of 
Figure  3)  If  it  is.  then  “F“  is  its  own  square  root;  therefore,  this  value  is 
returned  by  the  procedure,  (step  3  of  Figure  3).  If  "F"  was  neither  a  zero 
or  a  one.  then  it  is  checked  to  see  whether  it  is  negative.  Jf  it  is 
negative,  under  normal  circumstances,  an  error  would  be  assumed,  and  the 
procedure  would  terminate  rather  than  evaluating  for  a  complex  number. 
Due  to  the  nature  on  embedded  avionics  systems,  an  error  should  not  be 
fatal.  In  this  light,  rather  than  attempting  to  evaluate  a  complex  number, 
the  absolute  value  of  the  input  argument  is  formed,  (steps  A  and  5  of 


Figure  3)  The  built-in  function  for  absolute  value  was  found  to  give 
inconsistent  results,  so  the  absolute  value  is  found  by:  F  =  -  F.  If  this 
method  of  error  correction  proves  inappropriate  at  a  later  date,  it  would 
not  be  difficult  to  modify.  Perhaps  a  different  default  value  should  be 
assumed,  or  an  error  indicator  could  be  established. 

The  next  step  of  the  algorithm  is  to  obtain  the  exponent  portion  of 
the  input  argument,  (step  6  of  Figure  3)  JOVlAL's  specified  tables  make 
this  an  easy  conversion.  When  the  input  argument  was  placed  into  "F“,  "F" 
had  previously  been  established  as  a  table  whose  elements  identify  the 
components  of  the  floating-point  number.  Therefore,  the  item  “Feip"  is 
actually  the  exponent  portion  of  the  floating-point  number.  Immediately 
following  the  extraction  of  the  exponent,  this  same  exponent  portion  of 
the  floating-point  number  “F“  is  cleared  or  set  to  zero. 

The  next  few  steps  (steps  7  through  10  of  Figure  3)  are  to  compute  a 
polynomial  approximation  for  sqrt  (F).  Specifically,  the  computation 
begins  with  an  initial  approximation  of  #yQ"  shown  in  (22)  with 

successively  more  accurate  approximations  being  obtained  through  the  use 
of  Newton's  iteration  in  the  form  of  Heron's  formula. 


y0  -  .41731  +  .59016  *  f 


Yi  -  (Yi-i  *  f/Yi-i  )  /  2 
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The  coefficients  used  in  this  algorithm  are  those  presented  by  Cody 
and  Waite.  (4=  23)  The  approximation  described  by  Cody  and  Waite  is  in 
the  form  shown  in  (23).  Aside  from  the  original  calculation  of  "yq",  by 

examining  steps  8  through  10  of  Figure  3.  note  that  the  Newton  iteration 
is  performed  three  times.  Since  each  iteration  doubles  the  number  of 
correct  significant  digits  in  the  square  root,  this  assures  an  accuracy  of 
63.32  bits.  (4=  23)  The  next  step  is  to  determine  whether  or  not  the 
exponent  field  from  the  floating-point  number  originally  input  was  odd  or 
even,  (step  1 1  of  Figure  3)  Depending  upon  the  result  of  this  evaluation, 
different  actions  are  taken.  If  the  number  is  odd,  additional  calculations 
are  necessary  as  shown  in  equation  (21).  The  instruction  for  determining 
whether  the  number  is  odd  is  not  a  separate  function.  The  power  of  the 
I  JOVIAL  language  permit  testing  of  specific  bits.  Using  this  tool,  the 

low-order  bit  can  be  checked  to  determine  if  it  is  zero  or  one.  A  zero 
signifies  an  even  number  and  a  one  signifies  an  odd  number  which  is  just 
1  what  the  procedure  checks.  Given  that  multiplication  is  preferred 

multiplications  are  more  efficient  than  division  , 
the  calculation  sqrt(f)  /  sqrt(2)  is  represented  by  yj  *  sqrt(.5),  where 

the  decimal  representation  of  sqrt(.5)  is  the  constant  .707106781 1865. 

The  final  step  prior  to  returning  with  the  result  is  to  form  the 
exponent  portion  of  the  result,  (step  13  of  Figure  3) 
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There  are  three  steps  in  calculating  the  exponential  of  a 
floating-point  number.  The  first  step  is  the  reduction  of  the  given 
argument  to  a  related  argument  in  a  small  interval  symmetric  about  the 
origin.  The  second  step  is  the  computation  of  the  exponential  for  the 
reduced  argument,  and  the  final  step  is  the  reconstruction  of  the  desired 
function  from  its  components. 

The  exponential  is  formed  using  the  following  general  procedure.  Let 

X  -  N*ln(2)  ♦  g.  with  I  g  I  i  In  (2)/ 2  then 

eip(X)  *  eip(g)  *  CN  (24) 

The  accuracy  of  g  is  the  basis  for  the  accuracy  of  the  function  value. 
Let  y  -  eip  (g).  then  dy  /  y  -  dg.  This  means  that  the  relative  error  in 
eip  (g)  is  approximately  the  absolute  error  in  g.  This  error  is 
proportional  to  the  magnitude  of  X  when  X  is  exact  because  of  the 
finite  word  length  of  the  computer.  The  only  way  to  achieve  small 
absolute  error  in  g  is  to  extend  the  effective  precision  of  the  computer 
during  the  computation  of  g.  In  most  cases,  the  following  computation  is 


g  -  KXj  -  N*Cj)  +  X2l  -  N*C2  where 


Xj  +  x2  -x, 

X  i  is  the  integer  part  of  X  , 

Cj  +  C2  represents  In  (C)  to  more  than  working  precision 


This  method  gives  extra  digits  of  precision  equivalent  to  the  number 
of  extra  digits  in  the  representation  of  In  (0  when  N  is  small  enough 
that  N  *  Cj  is  representable  exactly  in  the  machine.  If  this  exact 

representation  cannot  be  accomplished,  the  computation  is  equivalent  to 

not  using  extra  precision.  Therefore,  the  magnitude  of  N  has  a  practical 

limit  which  results  in  a  limit  on  the  magnitude  of  X. 

* 

There  is  also  a  largest  and  smallest  X  such  the  exp(X)  can  be 
represented  in  the  machine.  For  example,  if  SMALLX  is  the  smallest 
positive  floating-point  number  and  BIGX  is  the  largest  without  causing 
overflow,  then  eip  (X)  can  be  represented  only  for  those  values  of  X 
that  between  In  (SMALLX)  and  In  (BIGX).  The  value  N  *  Cj  will  be 
representable  exactly  in  a  machine  for  any  X  within  the  specified  bounds, 
because  a  Cj  can  always  be  chosen  to  fit  the  bound.  Obviously,  careful 
argument  reduction  cannot  compensate  for  inaccuracies  in  X.  (4--  61) 


BIGX  88.02969193111 

SMALLX  -89.41598629223 

EPS  9.0949470 17729E- 13 

ONEOVERLN  1.4426950408890 

LN2  0.6931471805599 

Table  3  Constants  for  Exponential  Determination 

The  variable  “Arg“  is  the  argument  passed  to  the  square  root  function. 
Since  JOVIAL  treats  formal  arguments  as  read  only,  upon  program  entry, 
the  value  of  “Arg"  is  assigned  to  "X“.  (step  1  of  Figure  4)  The  variable 
“X"  is  then  used  throughout  the  remainder  of  the  procedure. 

The  constant  “BIGX“  (see  Table  3),  which  has  been  assigned  a  value 
that  is  slightly  less  than  the  natural  logarithm  of  the  largest  positive 
finite  floating-point  number  (step  2  of  Figure  4),  is  compared  with  the 
input  argument.  If  the  argument  is  larger  that  this  value,  an  error  would 
occur  during  calculating  its  exponential.  Since  this  application  is 
destined  for  embedded  avionics  systems,  a  solution  to  this  error  situation 
must  be  found  that  does  not  result  in  a  degradation  of  the  system.  The 
selected  solution  involves  replacing  the  input  argument  with  the  constant 
"BIGX“.  Obviously,  other  possible  options  are  available  to  resolve  the 
error  condition,  and  another  solution  can  easily  replace  the  existing 
methodology. 

The  constant  “SMALLX”  (see  Table  3),  which  has  been  assigned  a 
value  that  is  slightly  greater  than  the  natural  logarithm  of  the  smallest 
positive  finite  floating-point  number  (step  4  of  Figure  4),  is  compared 


Figure  4  Exponential  Structured  Flowchart 


with  the  input  argument.  If  the  argument  is  smaller  that  this  value,  an 
error  would  occur  during  calculating  its  exponential.  Again  the  discussion 
in  the  previous  paragraph  concerning  the  resolution  of  an  error  condition 
in  an  embedded  avionics  application  still  holds  true.  The  selected  solution 
involves  replacing  the  input  argument  with  the  constant  "SMALLX". 
Obviously,  other  possible  options  are  available  to  resolve  the  error 
condition,  and  another  solution  can  easily  replace  the  existing 
methodology. 

The  next  step  is  to  check  if  'X'  is  either  larger  than  a  positive  eps 
or  smaller  than  a  negative  eps.  (step  6  of  Figure  4)  In  either  case,  if  it 
is.  the  exponential  function  returns  a  value  of  1  and  terminates 
processing.  The  value  of  eps  (see  Table  3)  is  selected  with 
eip  (X)  =  J.O  to  machine  precision  such  that  |x|  <  eps  and  p,  *  X2 

will  not  underflow  for  |x|  i  eps.  Cody  and  Waite  have  suggested  that 
eps  =  2_t/2  where  there  are  t  base-2  digits  in  the  significand. 

The  next  step  (step  8  in  Figure  4)  involves  extracting  the  integer 
portion  of  the  floating-point  number  that  results  from  the  following  the 
calculation;  X  *  ll  /  in  (2)  ] .  As  noted  in  the  description  of  the  square 
root  function,  multiplication  is  not  as  costly  as  division.  Therefore,  the 
value  of  1  /  In  (2)  has  been  calculated  and  used  as  a  constant,  (see 
Table  3)  This  newly  formed  integer  is  then  transformed  into  a 
floating-point  number.  The  JOVIAL  specified  table  construct  is  put  to  use 
here  rather  than  calling  the  two  functions  INTRND  and  FLOAT.  Extraction 
of  a  specific  portion  of  a  floating-point  number  simply  involves  naming 
its  component  parts  and  using  these  names  to  access  the  needed  part. 


This  particular  construct  is  an  extremely  efficient  method  for  doing  this 
type  of  accessing,  and  is  not  confined  to  JOVIAL.  It  is  also  available  in 
the  Ada  language. 

The  computation  provided  by  Cody  and  Waite  that  is  specified  for  no 
guard  digits  is  use  to  create  a  new,  more  precise  number. 


g-l(  Xj  -  XN*C1)  +  X2]-XN*C2  where  (25) 

Xj  =  the  floating-point  value  of  the  integer  portion  of  X, 

x2=  x-xl§ 

q  -  0.693359375, 

C2  -  -2.121 9444005470E-4 

Now,  that  the  value  of  the  values  of  the  X's  and  N’s  are  known, 
equation  (25)  can  be  evaluated  for  g.  (step  9  in  Figure  4)  This  is 
followed  by  the  determination  of  the  rational  functions  R  (g)  which 
approximate  exp  (g)  /  2.  The  factor  of  0.5  is  inserted  to  counteract 
wobbling  precision.  The  calculation  of  the  coefficients  for  the 
approximation  are  determined  by  the  number  of  bits  in  the  significand. 
For  this  architecture,  the  number  of  bits  selected  are  between  30  and  42 
inclusively.  This  results  in  the  coefficient  list  of  Table  4  on  the  next 
page. 
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The  calculation  of  the  logarithm  required  three  steps.  First,  the 
given  argument  is  reduced  to  a  related  argument  in  a  small, 
logarithmically  symmetric  interval  about  one.  The  second  step  involves 
the  computation  of  the  logarithm  for  this  reduced  argument.  Finally,  the 
desired  logarithm  must  be  reconstructed  from  its  components. 
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Upon  entry  into  this  routine,  the  value  of  the  input  is  checked  to  see 
if  it  is  either  zero  or  less  than  zero,  (step  I  of  Figure  5)  If  it  is  either 
zero  or  negative,  under  normal  circumstances,  an  error  would  be  assumed, 
and  the  procedure  would  terminate.  In  this  function,  the  negative  value  of 
the  largest  floating-point  number  is  returned,  (step  2  of  Figure  5)  As 
previously  mentioned,  due  to  the  critical  nature  of  embedded  avionics 
systems,  an  error  should  not  be  fatal.  It  should  provide  an  alternate  path 
to  a  graceful  completion  of  the  function. 

Many  methods  exist  for  calculating  the  logarithm  of  a  reduced 
argument.  Cody  and  Waite  have  chosen  the  following  method.  (4,  42)  The 
initial  assumption  is  made  that  the  argument  is  in  the  following  form 

X  =  ±f  *  2e,  where  .5  s  f  <  1 

Determine  the  value  of  N  and  the  scaled  value  of  f  such  that 

X  -  f  *  2N.  where  .5  i  f  <  1 

initially,  f  is  assigned  the  value  of  the  input  argument.  This  allows 
for  modification  of  the  input  floating-point  number.  Then  an  estimate  for 
N  is  made.  N  is  given  the  value  of  the  exponent  of  the  input 
floating-point  number,  and  then  this  same  exponent  field  is  erased,  (step 
3  in  Figure  5) 

The  value  of  sqrt(.5)  has  been  previously  determined  and  stored  as  a 
constant  for  use  by  this  routine.  Depending  on  the  value  of  f ,  one  of  two 
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40 


aO  0.37339 1 68963 1 6E+ 1 

al  -0.6326086623386E+0 

a2  0.444455 1 5 1 0980E-2 

bO  0.4480700275574B+2 

bl  -0. 1 43 1 235435589E+2 

b2  0.1000000000000E+1 


Table  5  Coefficients  for  Polynomial  Approximation  to  Alog 


distinct  paths  may  be  taken.  The  value  of  f  is  compared  to  the 
sqrt(.5),  and  znum  and  zdem  will  vary  accordingly. 

After  forming  z  -  znum  /  zdem  and  w«z2,  evaluate  r(  z2)  -  w  * 
A(w)  /  B(w).  Both  A(w)  and  B(w)  are  polynomials  in  the  w  coefficients 
given  in  Table  5.  (step  8  in  Figure  5) 

Common  Logarithm  Implementation 

Obviously,  from  the  structure  chart  for  this  function  (Figure  6).  all 
the  work  is  done  by  the  Alog  function.  Since  the  JOVIAL  language  does 
not  support  multiple  entry  points,  the  common  logarithm  function  had  to 
be  formed  in  this  manner.  The  result  of  this  function  is  generated  through 
the  multiplication  of  the  natural  logarithm  of  the  input  argument  with  the 
natural  logarithm  of  “e“.  This  latter  item  is  encoded  as  a  constant  to 
avoid  wasted  effort  to  recalculate  for  every  use  of  this  function.  All  the 


Figure  6  Common  Logarithm  Structured  Flowchart 


restrictions  that  were  imposed  on  the  natural  logarithm  of  a  number  also 
apply  here. 


This  chapter  is  concerned  with  describing  the  methodology  used  for 
determining  the  correctness  and  performance  qualities  of  the  implemented 
functions.  Due  to  problems  in  the  availability  of  hardware  and  the 
associated  support  software,  the  testing  and  performance  evaluations  are 
somewhat  limited.  Hardware  became  available  towards  the  middle  of  the 
thesis  effort,  but  software  tools  used  for  development  were  incompatible 
with  those  required  by  the  available  1750A.  The  loader  used  by  the 
available  1750A  equipment,  expects  files  of  a  different  format  rhan  what 
is  created  by  the  software  development  tools.  Rather  than  developing  a 
new  loader,  a  routine  was  written  that  converts  load  modules  into  a 
format  required  by  the  1750A  loader.  The  reformatting  procedure  is  listed 
in  Appendix  C. 

Another  problem  that  had  to  be  overcome  before  testing  and 
evaluation  could  be  considered,  was  the  availability  of  input/output  (I/O) 
routines.  Without  I/O  routines,  further  considerations  for  testing  would  be 
fruitless.  No  I/O  packages  were  available,  and  as  a  consequence,  had  to  be 
created.  This  delayed  testing  efforts  considerably,  as  an  I/O  routine  had 
to  be  developed  with  the  use  of  the  N1L-STD-1750A  standard  ISA,  rather 


than  with  a  high-order  language.  The  1/0  package  developed  is  listed  in 
Appendix  B.  and  is  only  capable  of  writing  to  a  user  console. 

Performance  analysis  requires  the  comparison  of  1750A  results,  with 
those  generated  on  a  machine  of  higher  precision.  Unfortunately,  this 
requirement  made  the  newly  created  1/0  routine  insufficient  for  this  task. 
An  available  console  driver  has  a  routine  that  writes  user  specified  areas 
of  1750A  memory  to  magnetic  disk.  By  storing  a  function's  results  in  a 
specified  area  of  1750A  memory,  the  test  results  can  then  be  dumped  to 
disk  for  an  eventual  upload  to  a  VAX  11/780A.  The  results  are  then 
available  for  input  to  the  different  software  test  packages.  However,  the 
record  format  of  the  I750A  memory  dump  is  not  in  a  friendly  format,  and 
must  be  converted  to  a  readable  form.  At  the  time  of  this  writing,  a 
routine  for  making  the  disk  file  readable  is  not  completely  debugged. 
However,  it  is  at  a  point  where  it  could  be  completed  by  another 
programmer. 

The  aforementioned  problems  have  limited  the  amount  of  time 
available  for  designing  extensive  test  procedures.  Therefore,  validation, 
verification,  and  performance  analysis  is  confined  to=  manual  static 
analysis  methods,  critical  value  testing,  and  measurement  of  each 
algorithms  generated  error. 


Manual  Static  Analysis  Methods 


To  most  people,  manual  static  analysis  is  called  “desk  checking"  . 
Static  analysis  involves  the  search  for  any  inconsistencies  between  design 
tools  (i.e.  flowcharts),  design  details  (chapter  3),  program  headers,  and 
program  comments.  This  method  is  useful  for  finding  errors  caused  by  the 
translation  of  design  into  code,  as  well  as  possible  design  errors.  An 
inconsistency  may  indicate  potential  problems.  This  methodology  was 
used,  and  all  inconsistencies  that  were  found  were  resolved. 


Critical  Value  Testing 


Critical  value  testing  is  an  attempt  to  "break”  the  software,  and 
requires  the  selection  of  specific  arguments  that  could  possibly  cause 
problems.  A  knowledge  of  each  of  the  algorithms  is  required  to  select 
proper  arguments.  Individual  test  cases  are  not  listed  here,  but  the  reader 
may  find  specific  information  by  examining  the  test  procedures  listed  in 
Appendix  B. 

It  is  possible  to  generalize  the  tests  performed  without  listing  the 
specific  test  cases.  Potential  test  arguments  are  those  whose 
intermediate  results  could  generate  an  overflow  or  underflow,  or  are 
arguments  lying  in  the  fringe  of  computational  abnormality.  These 
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arguments  will  help  detect  problem  areas,  and  will  give  an  indication  as 
to  how  robust  each  function  is. 

In  addition,  arguments  that  test  each  path  of  the  algorithm  have  been 
selected.  Path  testing  is  limited  to  insuring  that  every  path  of  an 
algorithm  is  tested,  and  does  not  imply  that  every  possible  path 
combination  is  taken. 


Performance  Evaluation 

j  As  was  mentioned  in  the  introduction  of  this  chapter,  screen  output 

to  the  user  console  and  hard  copies  of  computed  results  are  insufficient 
for  performance  evaluation.  Their  use  would  imply  a  visual  comparison  of 
i  generated  results  against  published  tables.  Such  a  technique  limits  the 

number  of  comparisons  that  could  be  made,  and  would  cause  doubt  as  to 
the  credibility  of  the  comparisons.  At  best,  it  would  provide  a  good 
feeling  for  the  quality  of  each  function's  performance.  Therefore,  it  is 
better  to  automate  the  process  completely,  and  compare  the  generated 
results  against  another  machine  generated  standard. 

The  performance  evaluation  of  the  functions  involves  the  computation 
of  two  important  statistics:  the  maximum  relative  error  (MRE),  and  the 
root-mean-square  relative  error  (RE).  Their  values  are  determined  through 
the  use  of  (43)  and  (44),  where  F(x)  is  the  test  result  and  }(%)  is  the 
comparison  value  generated  by  the  same  extended-precision  function  call 
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written  for  the  VAX  1 1/780. 


This  method  of  error  checking  is  an  automatic  tabular  comparison, 
where  the  VAX  routines  serve  as  the  accepted  standard.  The  test  routine 
tests  densely  packed  samples  of  evenly  spaced  arguments  spread 
throughout  [-3tt,  3tt]  for  floating-point  algorithms,  and  [-1,1]  for 
fixed-point  algorithms.  When  regenerating  arguments  within  the  test 
modules,  it  is  important  not  to  introduce  unnecessary  errors.  This  means 
that  arguments  in  the  VAX  should  have  its  lower  order  bits  padded  with 
zeros.  The  most-significant  bits  must  be  equivalent  to  the  number  of  bits 
in  the  1750A  argument,  and  no  extra  precision  should  be  introduced. 

The  method  of  argument  generation  just  described  is  recommended  by 
Cody  (12;  762),  and  is  the  method  used  at  the  NASA  Lewis  Research 
Center.  This  method  is  preferred  to  a  random-number  test  because  it 
measures  the  relative  error  throughout  an  entire  interval.  Using  densely 
packed  arguments  also  gives  valuable  insight  to  problems  of  different 
argument  ranges,  if  the  evenly  spaced  interval  is  set  to  a  power  of  two 
(representable  on  both  machines),  and  is  not  less  than  the 


least-significant  bit  of  the  1750A  argument,  then  an  initial  argument  can 
be  chosen,  such  that,  zero  padding  will  only  have  to  be  performed  once. 
For  example,  if  an  initial  floating-point  argument  is  -3. 1415  and  the 

chosen  interval  is  2'2  .  the  second  argument  will  be  -3.1415  +  2‘2  . 

Additional  padding  is  not  necessary,  because  “carries“  are  cascaded 
forward  and  do  not  increase  the  number  of  most-significant  bits  in  the 
next  argument.  Arguments  used  in  the  function  calls  on  both  machines 
must  be  the  same,  and  must  be  generated  in  the  same  order. 

Extra  care  is  needed  while  reading  the  I750A  results  from  disk.  Each 
of  the  1750A  results  are  stored  in  an  unformatted  file,  and  must  be  read 
into  a  binary  record.  This  record  is  moved,  bit-by-bit,  to  a  variable  of  the 
appropriate  type  (VAX  11/780  fixed-point  or  floating-point).  The 

bit-by-bit  manipulation  is  accomplished  through  the  use  of  JOVIAL 
specified  tables,  and  prevents  conversion  errors  associated  with 

formatted  input. 

Before  a  comparison  of  the  two  results  (one  from  the  1750A,  and  the 
other  from  the  VAX)  can  be  made,  the  results  generated  within  the  test 
module  must  be  reduced  to  the  same  precision  (same  number  of 
most-significant  bits)  as  those  from  the  1750A.  The  precision  reduction 
gives  a  rounded  result  that  can  be  used  to  determine  the  MREand  RE,  and 
will  give  a  meaningful  interpretation  to  the  inherited  error  of  the  1750A 
functions. 
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The  purpose  of  this  thesis  was  to  develop  and  to  do  performance 
evaluation  on  a  run-time  math  library  developed  specifically  for 
M1L-STD-1750A  architectures.  The  library  consists  of  the  floating-point 
implementation  of  several  algebraic  functions.  Performance  evaluation 
was  the  major  effort  of  this  thesis,  but  not  in  the  manner  intended. 

Function  approximations  are  accomplished  through  the  use  of  either 
Chebyshef  or  rational  approximations.  The  two  different  approximation 
methods  were  discussed  in  chapter  two,  and  are  useful  in  understanding 
certain  design  considerations.  The  values  of  each  polynomial's 
coefficients  were  derived  by  (or  were  modifications  of  those  derived  by) 
Cody  and  Waite.  (4;  17-84)  However,  the  implementation  designs  are 
significantly  different  from  those  suggested  by  Cody  and  Waite.  The 
primary  difference  between  the  implemented  designs  and  those  suggested 
by  Cody  and  Waite,  are  the  methods  of  argument  reduction  required  of  each 
function. 

Performance  evaluation  turned  out  to  be  the  major  effort,  but  not 
because  of  extensive  or  elaborate  testing  of  the  library  functions.  Most  of 
the  effort  involved  overcoming  the  following  problems; 


1. )  There  were  several  compiler  bugs  in  the  original  1750A  compiler 
used.  Assembly  listings  had  to  be  reviewed,  in  order  to  verify  each 
compilation  of  the  source  code. 

2. )  The  use  of  a  simulator  for  performance  evaluation  was  ruled  out 
because  of  the  limited  number  instructions  that  could  be  simulated,  its 
inability  to  simulate  the  use  of  floating-point  data,  and  the  relative  speed 
at  which  results  were  calculated.  The  simulator  also  lacked  a  facility  for 
writing  results  to  mass  storage.  Storage  of  results  on  an  external  device 
is  necessary  for  input  to  software  test  packages. 

3. )  A  new  compiler  and  linker  was  introduced  near  the  midpoint  of  the 
thesis  effort,  and  required  a  long  learning  curve  in  order  to  use  them. 

4. )  Once  a  I750A  machine  became  available,  it  was  determined  that 
all  its  support  software  was  intended  for  use  with  files  created  by  the 
old  compiler  and  linker. 

5. )  Rather  than  use  a  compiler  and  linker  that  had  several  deficiencies, 
or  write  a  new  loader  routine,  it  was  decided  to  write  a  support  tool  that 
would  convert  load  modules  into  a  format  expected  by  the  available  loader. 

6. )  The  reformatting  program  required  the  use  of  JOVIAL  and  its 
specified  table  features.  It  also  required  the  use  of  FORTRAN  routines  to 
perform  the  I/O  of  source  and  target  files.  The  FORTRAN  and  JOVIAL 
interfaces  did  not  operate  as  expected,  and  the  use  of  COnnON/COMPOOL 
areas  wouldn't  work.  This  required  parameter  passing  between  the 
routines,  and  the  documentation  for  this  type  of  interface  was  very 
inadequate;  however,  the  problems  were  eventually  resolved. 


7. )  The  reformatting  tool  was  written  for  use  on  a  VAX  1 1/780.  It 
was  assumed  that  the  JOVIAL  compiler  was  free  of  bugs  for  a  VAX  target. 
However,  when  the  reformat  routine  was  being  debugged,  it  was 
discovered  that  JOVIAL  table  names  could  be  overlayed,  but  corresponding 
table  items  weren’t  overlayed  with  them.  This  problem  took  a  long  time 
to  discover,  and  an  additional  amount  of  time  to  design  around. 

8. )  I/O  routines  have  not  been  written  for  the  1750A,  and  had  to  be 
developed.  These  routines  are  only  capable  of  writing  to  a  console  screen. 

9. )  Screen  output  is  insufficient  for  generating  the  thousands  of 
results  that  would  be  needed  during  testing  and  evaluation,  so  another 
means  of  capturing  the  data  had  to  be  developed.  Due  to  the  lack  of  time 
and  inexperience  in  the  internal  I/O  communications  techniques  of  the 
I750A  hardware,  development  of  a  disk  I/O  routine  was  not  a  feasible 
alternative.  It  was  determined  that  results  could  be  stored  in  specific 
locations  of  memory,  and  then  an  available  console  routine  could  be  used 
to  write  the  information  to  disk.  An  additional  problem  was  encountered 
when  it  was  discovered  that  the  record  format  of  the  disk  files  is  not  in  a 
VAX  friendly  format,  and  another  routine  had  to  be  written  to  unpack  the 
stored  results. 

These  problems  limited  the  scope  of  this  thesis  effort  to  developing 
the  following:  designs;  code  that  is  free  of  syntax  errors;  the 
development  of  command  files  for  compiling,  assembling,  and  linking 
routines  written  for  the  1750A;  tools  for  formatting  load  modules  that 
are  capable  of  being  loaded  into  a  Sperry  1631  implementation  of  the 
M IL-STD- 1 750A;  and  tools  that  unpack  test  results  stored  on  an  RT/11 
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formatted  floppy  disk.  Generic  test  algorithms  are  provided,  but  are  not 
written  in  a  high-order-language.  They  provide  the  basic  structure  for 
critical  range  testing,  and  a  means  of  evaluating  and  measuring  each 
functions  performance. 


Recommendations 


The  products  produced  by  this  thesis  effort  are  at  point  where  design 
of  the  intended  performance  evaluation  can  begin.  All  the  groundwork  has 
been  provided,  and  should  be  adequate  for  someone  to  continue  the  effort. 
Many  of  the  aforementioned  problems  have  been  resolved,  and  support  tools 
and  command  files  are  provided  to  shorten  the  learning  curve  that 
follow-on  programmers  will  have  to  experience. 

The  following  recommendations  should  be  considered  if  this  effort  is 
continued. 

1. )  If  the  effort  is  limited  to  the  use  of  JOVIAL,  an  analysis  should  be 
made  for  determining  how  to  handle  exceptions  detected  at  run  time. 
Exceptions  include  arguments  outside  legally  defined  limits. 

2. )  Since  Ada  has  features  for  exception  handling,  all  the  library 
functions  should  also  be  developed  and  implemented  in  Ada. 

3. )  Another  point  may  be  in  favor  or  using  Ada  is  that  it  also  allows 
the  creation  of  generic  packages  and  subprograms.  The  generic 
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subprograms  define  a  template,  and  generic  parameters  provide  the 
facility  for  tailoring  the  template  to  fit  a  particular  need  at  translation 
time.  In  other  words,  one  subprogram  could  provide  calculations  for  both 
fixed-point  or  floating-point  arguments,  based  on  how  it  is  used  at 
compile  time.  Because  a  generic  package  would  not  be  able  to  take 
advantage  of  the  specific  hardware  functions  unique  to  floating-point  and 
fixed  point  routines,  this  may  result  in  a  degradation  of  performance. 

4. )  Initially,  it  was  discussed  that  all  the  math  library  routines  should 
be  written  in  both  JOVIAL  and  Ada  with  the  intent  that  a  comparative 
evaluation  could  be  done  on  the  two  languages.  Unfortunately,  an  Ada 
compiler  targeted  to  the  1750A  is  not  yet  available.  When  a  compiler 
does  become  available,  it  is  recommended  that  a  new  Ada  math  library  be 
developed  and  this  comparative  evaluation  be  performed. 

5. )  The  compiler  problems,  mentioned  above,  should  be  corrected,  and 
I750A  architectures  and  associated  support  software  should  be  acquired 
before  more  time  is  allocated  to  the  effort. 


7^.V.V  .’W*  W 


Appendix  A 


The  following  pseudo-operations  were  used  in  describing  the 
implementation  designs  of  the  different  mathematic  functions. 

ADX(X.N):  augments  the  integer  exponent  of  a  floating-point 

representation  of  X  by  N.  This  scales  the  argument  X  by  2N. 
For  example. 

ADX(  1.0.2)  =  4.0 

FIX(X):  returns  the  fixed-point  representation  of  the 

floating-point  value  X  .  This  operation  requires  explicit 
conversion  in  JOVIAL. 

FLOAT(X):  returns  the  floating-point  representation  of  the 

fixed-point  argument  X.  This  operation  requires  explicit 
conversion  in  JOVIAL. 


ODD(X):  determines  whether  the  argument  X  is  odd.  For  an 

integer,  the  least-significant  bit  is  checked  directly.  For  a 
floating-point  number,  the  integer  portion  is  checked.  A 
description  of  the  floating-point  process  for  this  determination 
is  given  below. 


Figure  7  Bit  Layout  of  1750A  Floating-Point  Number 


To  determine  whether  the  integer  portion  is  odd,  knowledge 
of  the  internal  representation  of  the  1750A  floating-point  number  is 
necessary.  The  argument  X  is  a  JOVIAL  specified  table  item  that 
makes  the  components  shown  in  Figure  7  easily  accessible.  Within 
this  table  is  an  integer  item  that  overlays  the  exponent  field  of  X. 
This  exponent  field  is  the  tool  needed  to  check  whether  the  integer 
portion  is  odd  or  even.  Since  X  has  a  value  of  one  or  greater,  and 
all  floating-point  values  are  normalized,  the  exponent  can  be  used  to 
point  to  the  least  significant  bit  of  the  integer  field.  Because  X  is 
positive,  a  one  in  the  least  significant  bit  would  indicate  the 
integer  portion  is  odd.  A  limit  on  the  maximum  value  of  the 
coefficient  has  been  imposed  by  the  functions  that  use  this  routine. 
This  limit  prevents  the  least-significant  bit  of  the  integer  portion 
from  falling  in  the  exponent  or  "LSB”  area  of  the  floating-point 
coefficient  (see  Figure  7). 


Since  the  1750A  architecture  requires  that  all  floating-point 
values  be  normalized,  the  most-significant  bit  is  in  the  first  bit 
position  following  the  sign  bit.  The  decimal-point  is  assumed  to  be 
positioned  immediately  behind  the  sign  bit,  but  immediately  in  front 
of  the  most-significant  bit.  The  exponent  represents  a  power  of 
two;  therefore,  if  e  represents  the  value  of  the  exponent  field,  the 
value  of  the  floating-point  number  is:  coefficient  *  2e  . 
Equivalently,  it  is  obvious  that  the  decimal-point  floats  e  places  to 
the  left  if  e  is  negative,  or  s  places  to  the  right  if  positive. 

Knowledge  of  how  floating-point  numbers  are  stored  can  be 
used  to  determine  whether  the  integer  portion  of  a  number  is  odd. 
The  following  example  gives  an  explanation  of  the  process. 

Given  the  following  machine  representation  of  a 

% 

floating-point  number,  determine  whether  its  integer  portion 
is  odd.  In  the  example  below,  the  decimal-point  was  inserted 
only  for  clarity. 

0. 1 1 0000000000000000000000000000 1 0000000000000000 

Since  the  sign  bit  of  the  exponent  is  zero,  the  value  of 
the  coefficient  is  positive.  The  following  two  numbers  are 
summed  together  to  determine  the  value  represented  by  this 
coefficient: 


The  exponent  field  is  in  bold  text,  and  has  the  value 
one.  Therefore,  the  value  of  this  floating-point 
representation  is,  the  coefficient  (.75)  multiplied  by  two 
to-the-power-of  the  exponent  ( 1 ),  or  1 .5. 

.75  *  21  =  1.5 


Another  way  to  compute  the  result  is  to  shift  the 
decimal-point  in  a  direction  as  indicated  by  the  exponent. 
The  exponent  in  this  case  is  +1  ,  so  the  decimal-point  is 
shifted  one  position  to  the  right.  The  number  can  then  be 
computed  in  a  similar  manner  as  described  above. 

This  last  method  demonstrates  how  to  determine 
whether  this  example  is  even  or  odd.  If  the  decimal-point 
is  shifted  I  position  to  the  right,  this  number  will  have  I 
integer  bit  and  38  fractional  bits.  The  integer  bits  always 
occupy  the  left-most  position  of  the  number.  If  the  exponent 
is  thought  of  as  a  pointer  from  the  left-most  side  of  the 
number,  the  least-significant  integer  bit  can  be  found.  The 
exponent  in  this  example  points  to  bit  position  one.  Since 
the  bit  is  set  to  1,  this  example's  integer  value  is  odd.  I 


■L  I  "Ji 


iii|i  in  wi  i  '»j  ■  w.1  \m  ^ppr^  V*  r^y^ywT 


INT(X):  return  the  integer  portion  of  the  floating-point 

argument  X.  The  description  ODD(X)  given  above  determines  the 
least-significant  bit  of  the  integer  portion  of  the  floating-point 
argument.  This  is  used  to  extract  the  entire  integer  portion  of 
the  argument  (bits  0  through  the  teast  significant  bit). 
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DATE: 

VERSION: 

MfiMC  . 

rVTE: 

HOOULE  NUMBER: 
DESCRIPTION: 


30  August  1085 
1.0 
ALog 
1.0 


This  function  Is  cal  ltd  to  compute  the  natural  log  of 
of  tha  arguMnt  'Arg‘.  Sine* 

ALog10<X>  -  ALog«>  *  ln<e> 

ALog  is  also  cat  lad  bg  ALog 10  to  do  its  computations. 
PASSED  VARIABLES:  Arg  -  an  axtandod  precision  floating-point  variable 


RETURNS: 
NOOULES  CALLED: 
AUTHOR: 

HISTOAY: 


Tha  natural  log  of  arg  in  axtandod  precision  float 


Copt.  Jennifer  Fried 

This  project  aas  undertaken  as  a  thesis  project  for 
partial  fulfillment  of  requirements  for  an  NS  degree 
in  Information  Science  from  the  Air  Force  Institute 
of  Technology.  Sponsoring  organization  is  the  ASD 
Language  Control  Branch,  Wright  Patterson  AFB,0h. 


***> 


STAAT 

OEF  PAOC  ALog  RENT <  Arg  )  F  30; 
BEGIN 


ITEH 

Arg 

F 

30; 

ITEM 

Nn 

S 

7; 

ITEH 

Xn 

F 

39; 

ITEM 

Znum 

F 

39; 

ITEH 

Zden 

F 

39; 

ITEH 

Zz 

F 

39; 

ITEH 

Rz 

F 

39; 

ITEH 

Az2 

F 

39; 

ITEH 

We 

F 

39; 

ITEH 

Ae 

F 

39; 

ITEH 

Be 

F 

39; 

TABLE  Over lays  <0>  U  3; 


BEGIN 

ITEH 

Ff 

F  39 

P0S(0,0>: 

ITEH 

Fexp 

S  ? 

P0S<8, 

END 

CONSTANT 

ITEH 

Zero 

F 

39 

m 

0.0; 

CONSTANT 

ITEH 

PtFive 

F 

30 

m 

0.5; 

CONSTANT 

ITEH 

SqrtPtFive  F 

39 

m 

0.7071067811865; 

CONSTANT 

ITEH 

AO 

F 

39 

m 

0.37339 168963 18E+1; 

CONSTANT 

ITEH 

A1 

F 

39 

m 

-0.6326086623386E+0; 

CONSTANT 

ITEH 

A2 

F 

39 

M 

0.444455 13 10980E-2; 

CONSTANT 

ITEH 

BO 

F 

39 

m 

0 . 4480700275574E+2 ; 

CONSTANT 

ITEH 

B1 

F 

39 

m 

-0.  I43I235433589E+2; 

CONSTANT 

ITEH 

B2 

F 

39 

m 

0.  1000000000000E+1; 
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OflTE: 

30  August  1985 

UERSION: 

1.0 

IWc! 

RLog  10 

N00ULE  NUMBER: 

1.0 

DESCRIPTION: 

This  subroutine  Is  called  to  coepute  the  base  10  log 

of  the  passed  argueent.  Since 

ALoglO  ■  RLog  *  log(e) 

It  cakes  a  call  to  RLog 

PASSED  VARIABLES:  firg  -  an  extended-precision  floating-point  variable 

RETURNS: 

The  floating-point  rep  of  log(Arg) 

NOOULES  CALLED 

RLog 

AUTHOR: 

Capt.  Jennifer  Fried 

HISTORY: 

This  project  eas  undertaken  as  a  thesis  project  far 

partial  fulfil  leant  of  requl resents  for  an  NS  degree 

in  Information  Science  from  the  Air  Force  Institute 

of  Technology.  Sponsoring  organization  is  the  ASO 

Language  Control  Branch,  Wright  Patterson  AFB,0h. 

START 


REF  PROC  RLog  RENT  <  firg  )  F  39; 

BEGIN 

ITEM  Arg  F  39; 

END 

DEF  PROC  RLog 10  RENT  <flrg>  F  39; 

BEGIN 

ITEM  ft»g  F  30; 

CONSTANT  ITEM  Log'*'  F  30  -  0.4342944819033; 


END 

TERN 


RLog 10  ■  ALog(Arg)  *  Log'*'; 
RETURN; 
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DATE;  30  August  1985 

VERSION:  1.0 

NAME:  Exp 

MOOULE  HUMBER:  1.0 
DESCRIPTION 

Returns  the  extended-precision  floating-point  value 
for  e**Arg 

PftSSED  URRIflBLES:  Rrg  -  an  extended  precision  floating-point  variable 

RETURNS:  e**x 

HOOULES  CflLLED:  none 

AUTHOR:  Copt.  Jennifer  Fried 

HISTORY:  This  project  eas  undertaken  as  a  thesis  project  for 

partial  fulfil leent  of  requireeents  f or  an  HS  degree 
in  Inforeation  Science  froe  the  Air  Force  Institute 
of  Technology.  Sponsoring  organization  is  the  RSD 
Language  Control  Branch,  Uright  Patterson  HFB,0h. 


START 

OEF  PAOC  Exp  RENT  <Arg>  F  39; 
BEGIN 

ITEM  RRO  F  39; 


ITEM  Xx 

F 

39; 

ITEM  Xn 

F 

39; 

ITEM  Og 

F 

39; 

ITEM  XI 

F 

39; 

ITEM  X2 

F 

39; 

ITEM  Zz 

F 

39; 

ITEM  Pz 

F 

39; 

ITEM  Qz 

F 

39; 

ITEM  Nn 

S 

7; 

TABLE  Overlays 

<0>  U  3 

i 

BEGIN 

ITEM 

Bg 

F  39 

P0S(0,0) 

ITEM 

Rexp 

S  7 

P08<8, 1 ) 

END 

CONSTANT 

ITEM 

Xoax 

F 

39  - 

1.701411834599E+38; 

CONSTANT 

ITEM 

Xein 

F 

39  > 

1.469367938527E-39; 

CONSTANT 

ITEM 

Xbig 

F 

39  > 

88.02969193111; 

CONSTANT 

ITEM 

Xseall 

F 

39  - 

-89.41598629223; 

CONSTANT 

ITEM 

Eps 

F 

39  - 

9. 0949470 17729E- 13; 

CONSTANT 

ITEM 

PO 

F 

39  - 

0 . 2499999999999E+0; 

CONSTANT 

ITEM 

PI 

F 

39  - 

0.5950425497759E-2; 

CONSTANT 

ITEM 

00 

F 

39  - 

0 .50000000000006*0; 

CONSTANT 

ITEM 

Q1 

F 

39  - 

o! 535675 1764522E-1; 

CONSTANT 

ITEM 

02 

F 

39  - 

0 . 2972936368224E-3 ; 

CONSTANT  ITEM  Cl 
COMSTflNT  ITEM  C2 


F  39  *  0.693359375; 

F  39  ■  -2. 1219444005470E-4 


COMSTflMT  ITEM  Ln2  F  39  ■  0.0931471805599; 

COMSTflMT  ITEM  0n*0w*rLn2  F  39  -  1.4426950408890; 

COMSTflMT  ITEM  On*  F  39-1.0; 

COMSTflMT  ITEM  PtFlv*  F  39  ■  0.5; 


Xx  -  Arg; 

IF  f*g  >  Xbig; 

Xx  ■  Xbig; 

IF  Xx  <  Xsaall; 

Xx  ■  Xsaall; 

IF  (Xx  <  Eps)  AND  (Xx  >  -Eps); 

Exp  ■  On*; 

ELSE 

BEGIN 

Mn  ■  (*  S,  fl  7  *><  Xx  *  Qn*0u*rLn2  >; 

Xn  ■  (*  F  39  *)(  Nn  ); 

XI  ■  (*  F  39  *)((*  S  31  *)(  Xx  »; 

X2  *  Xx  -  XI; 

Og  »  ((XI  -  Xn  *  Cl)  ♦  X2>  -  Xn  *  C2; 

Zz  ■  Og  *  Og; 

Pz  ■  (PI  *  Zz  ♦  PO)  *  Gg; 

Qz  ■  (Q2  *  Zz  ♦  Q1>  *  Zz  ♦  QO; 

Rg(0>  -  PtFiv*  ♦  Og  *  Pz  /  (Qz  -  Pz>; 

R*xp(0>  ■  R*xp(0>  +  Nn  ♦  (*  S  7  *>(  1  >; 

Exp  ■  flg(0>; 

END 

HETUflM; 
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OflTE:  30  August  1065 

UER8I0N:  1.0 

NRHE:  Sqrt 

MOOULE  NUMBER:  1.0 
DESCRIPTION: 

Approximates  the  square  root  of  th  argues nt  'Ary* 
PASSED  URRIRBLES:  Rrg  -  an  extended  precision  floating-point  uariefch 


RETURNS: 

NODULES  CALLED: 
AUTHOR: 

HISTORV: 


An  extended  precision  float  representation  of  the  sqrt 
of  ‘Arg’ 


Copt.  Jennifer  Fried 

This  project  eas  undertaken  as  a  thesis  project  for 
partial  fulfil  leant  of  requ I resents  for  an  NS  degree 
in  Information  Science  froa  the  Air  Force  Institute 
of  Technology.  Sponsoring  organization  Is  the  RSO 
Language  Control  Branch,  Wright  Patterson  AFB,0h. 


ORTE:  19  JULY  1969 

UER6I0N:  1.0 

NAME:  llathLib 

MODULE  NUMBER:  1.0 
DESCRIPTION: 

This  cospool  is  rsquirsd  by  any  JOUIRL  progras  that  nssds  to 
rsfsrsncs  any  of  the  sath  functions  srittan  for  floating-point 
or  fixed-point  cosputations 
PRSSED  UARIABLES:  N/R 


RETURNS: 
MODULES  CALLED: 
AUTHOR: 


HISTORY: 


Capt.  Steven  R.  Hotchkiss  and 
Ctpt  Jsnnifsr  Frisd 

This  project  sas  undertaken  as  a  thesis  project  for 
partial  fulfil  leant  of  requireeents  for  an  NS  degree 
in  Information  Science  from  the  Rir  Force  Institute 
of  Technology.  Sponsoring  organization  is  the  ASD 
Language  Control  Branch,  Uright  Patterson  AFB,0h. 


-  •  v*‘  «  I 


wA 


-*•  ..*•  -V 


•  A'  .j 


PflOC  Tan  RENTCXx)  R  12,18 

BEGIN 

ITEM  Xx  R  1,30; 

END 

PflOC  Cot  RENTCXx)  R  12,18 

BEGIN 

ITEM  Xx  R  1,30; 

END 

PflOC  RSin  RENTCXx)  R  1,30 

BEGIN 

ITEN  Xx  fl  1,30; 

EM) 

PflOC  flCos  RENTCXx)  R  1,30 

BEGIN 

ITEN  Xx  fl  1,30; 

END 

PflOC  RTan  RENTCXx)  R  1,30 

BEGIN 

ITEN  Xx  fl  1,30; 

END 

PflOC  Sinf  RENTCXx)  F  38 

BEGIN 

ITEN  Xx  F  38; 

END 

PflOC  Cosf  RENTCXx)  P  38 

BEGIN 

ITEN  Xx  F  38; 

END 

PflOC  Tanf  RENTCXx)  F  38 

BEGIN 

ITEN  Xx  F  38; 

END 

PflOC  Cotf  RENTCXx)  F  38 

BEGIN 

ITEN  Xx  F  30; 

END 

PflOC  ftSInf  RENTCXx)  F  38 

BEGIN 

ITEN  Xx  F  38; 

ENO 

PflOC  RCotf  RENTCXx)  F 

BEGIN 

ITEN  Xx  F  30; 
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OflTE: 

UERSION: 

NAME: 

nOOULE  HUMBER: 
DESCRIPTION: 


20  August  1085 

1.0 

loRsfs 

1.0 


This  Coapool  Is  necessary  to  reference  routlnss  that 
necessary  for  tasting  and  psrforaanea  evaluation  of  ail  aath 
functions  developed  for  tha  1750. 

PASSED  VARIABLES:  N/A 


RETURNS: 

MOOULES  CALLED: 
AUTHOR: 

HISTORY: 


N/A 

H/H 


and 


Copt.  Steven  A.  Hotchkiss 
Copt.  Jarri far  Friad 

This  projact  sos  undertaken  as  a  thesis  project  far 
partial  fulfil laant  of  requirements  far  an  NS  degree 
in  Information  Science  fraa  tha  Air  Force  Institute 


of  Technology.  Sponsoring  organization  is  the  ASD 
Language  Control  Branch,  Height  Patterson  AFB,0h. 


AT 


CONPOOL  loAefs; 


*  The  fol losing  ITEM*  are  required  to  print  a  carriage  return  and 
'  line  feed  on  a  terminal  connected  to  a  HIL-STD-1750  computer 

m 


DEF  ITEM  Carriage  STATIC  U  16  -  2573; 

OEF  I  TEH  CRLF  STATIC  C  2; 

OUERLAV  Carriage:  CALF; 


The  fol losing  referenced  subroutine  is  sritten  in  1750  Assembly  language 
and  is  used  to  print  character  strings  only.  Honcharactar  types  sill 
have  to  be  converted  before  calling  this  routine.  The  fol losing  DEFINE  is 
recoeaended  for  all  routines  calling  ObcSia: 

DEFINE  URITE'STRINO(R>  '  ,Printc(UOROSIZE<IR),LOC(!R»"; 

An  example  of  a  typical  call  foiloss: 

ITEM  Example  C  2; 


UR  I TE ' STfi I N0<Examp I e  >; 


PROC  Printc  RENT (Length,  Message); 


ITEM  Langth  U  (BITSIMU0RD-1); 
ITEM  Massaga  P; 

END 


Tha  fol lowing  rafarancad  routina  is  nacassorg  for  routinas  aishing 
to  conuart  floating-point  uaiuas  to  a  character  string 


■PROC  FltToChar  (flrg)  C  20; 
BEGIN 

ITEM  Arg  F  39; 

END 


Tha  fol losing  rafarancad  routina  is  nacassorg  for  routinas  aishing 
to  congart  fixad-point  galuas  to  a  charoctar  string.  Tha  gcriabla 
IntOgarlag  aust  ba  ouarlagad  on  top  of  a  fixad-point  gcriabla  and 
BitsInFrac  is  an  intagar  galua  indicating  tha  nuabar  of  fractional 
bits  in  tha  fixad-point  galua. 


EF  PROC  FlxToChcr  (IntOgarlag,  BitsInFrac)  C  20; 
BEGIN 

ITEM  IntOgarlag  8  31; 

ITEM  BitsInFrac  U  8; 


DATE:  29  August  1985 

UERSION:  1.0 

NAME:  FixToChcr 

MODULE  HUMBER:  1.0 

DESCRIPTION: 

This  rout  I ns  is  ussd  to  convert  fixed-point  values  into 
character  representation.  This  rout I  ns  sets  necessary  for 
testing  and  performance  svaluation  of  math  routines  developed 
for  the  1750 

PASSED  UARIABLES:  lnt0v«May  -  An  Integer  Variable  Overlayed  on  top 

of  a  f ixed-ooint  value 

BltsinFrac  -  the  nueber  of  fractional  bits  of 
the  fixed-point  argument 

RETURNS:  a  20  character  representation  of  the  argument 

MODULES  CALLED:  FltToChar 

AUTHOR:  Copt.  Steven  A.  Hotchkiss  and 

_  Capt.  Jennifer  Fried 

HISTORV:  This  project  eas  undertaken  as  a  thesis  project  for 

partial  fulfil  leant  of  requirements  for  an  MS  degree 
in  Information  Science  from  the  Air  Force  Institute 
of  Technology.  Sponsoring  organization  is  the  ASD 
Language  Control  Branch,  Urlght  Patterson  RFB,0h. 


RETURNS: 
MODULES  CALLED: 
AUTHOR: 

HISTORV: 


3fc  I 


FixToChar  -  FltToChar<Arg<0» 
RETURN: 


DOTE:  29  August  1985 

UERSIQN:  1.0 

NAME:  FltToChor 

MODULE  NUMBER:  1.0 

DESCRIPTION: 

This  rout ins  is  usad  to  convert  floating-point  values  into 
charoctar  raprasantation.  This  routina  aas  necessary  for 
tasting  and  perforeance  evaluation  of  oath  routinas  developed 


for  tha  1750 

PASSED  URRIABLES:  Arg  -  tha  value  to  ba  converted 

RETURNS:  a  20  charoctar  raprasantation  of  tha  orguaant 

MOOULES  CALLED:  nona 

AUTHOR:  Copt.  Stavan  A.  Hotchkiss  and 

Copt.  Jannl far  Fr lad 

H I STORY :  This  pro j act  aas  undartakan  as  a  thasis  project  for 

partial  fulfil laant  of  raquireaants  for  an  MS  dagraa 
in  Inforaation  Scianca  froa  tha  Air  Forca  Institute 


of  Technology.  Sponsoring  organization  is  the  ASO 
Language  Control  Branch,  Wright  Patterson  AFB,0h. 


START 


DEF  PROC  FltToChor  (Arg)  C  20; 
BEDIM 


DEFINE  Vas 

"18*1 

DEFINE  No 

-1B*0 

ITEM  ftrg 

F 

39; 

ITEM  Fraction 

F 

39; 

ITEM  Teap 

F 

39; 

ITEM  Result 

C 

20; 

ITEM  lx  U  8; 

ITEM  ly  U  8; 

ITEM  ExpCnt  U  8; 

ITEM  NagExp  8; 

ITEM  CharUal  U  8; 

ITEH  CharAep  C  1; 

OUERLAV  CharAep:  CharUal; 

ITEM  ZoroRep  STATIC  C  1  -  'O'; 

ITEM  ZaroUal  STATIC  U  8; 

OUERLAV  ZoroRep:  ZaroUal; 

CONSTANT  ITEM  Zero  F  39  ■  0.0; 

CONSTANT  ITEM  One  F  39  -  1.0; 

CONSTANT  ITEM  TanFloat  F  39  -  10.0; 

CONSTANT  ITEM  PtFive  F  39  -  0.3; 

CONSTANT  ITEM  PtOna  F  39-  0.1; 
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Rasul t  -  '  O.OOOOOOOOOOOOOE+OO'; 

IF  flrg  <  Zaro; 

BEGIN 

Fraction  ■  -flrg; 

BYTE (Rasul t,0, 1  >  - 
END 
ELSE 

Fraction  ■  flrg; 

IF  Fraction  <  PtOna; 

NagExp  ■  Vas; 

ELSE 

NagExp  ■  No; 

ExpCnt  -  0; 

WHILE  (Fraction  >  Ona); 

BEGIN 

ExpCnt  »  ExpCnt  +  I; 

Fraction  ■  Fraction  /  TanFloat; 

END 

IF  (NagExp  ■  Vas)  AND  <  Fraction  <>  Zaro); 

BEGIN 

BVTE(Rasult, 17,1)  • 

WHILE  (Fraction  <  PtOna); 

.  BEGIN 

%  0  ExpCnt  «  ExpCnt  +  1; 

Fraction  «  Fraction  *  TanFloat; 

END 

END 

iy  ■  0; 

WHILE  ((Fraction  <>  Zaro)  AND  (ly  <  13>>; 

BEGIN 

Taap  ■  Fraction  *  TanFloat; 

IF  ly  ■  12; 

Taap  ■  Taap  +  PtF iva; 

Chortle  I  -  (*  U  8  ♦)(  Taap  >; 

Fraction  ■  Taap  -  (*  F  30  *>(  Chordal  ); 

Chorda I  ■  Chorda I  +  Zarodal; 

BYTE (Rasul t, ly+3, 1 )  ■  Charflap; 

iy  ■  iy  +  t; 

END 

Chordal  ■  (*  U  8  *)(ExpCnt  NOO  10)  +  Zarodal 

BYTE (Rasul t, 19, 1)  -  Charflap; 

Chordal  ■  (*  U  8  ♦XExpCnt  /  10)  +  Zarodal 

BYTE(Rasult, 18, 1)  -  Charflap; 

FltToChar  ■  Rasul t; 

RETURN; 

END 

TERN 
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TITLE  HQL<PRINTC> 

nOOULE  PRINTC 


DOTE:  4  Septeeber  1989 

VERSION:  1.0 

NAME:  Printc 

nOOULE  MUHBER:  1.0 
DESCRIPTION: 

This  eodule  is  called  to  print  a  character  string  onto 
a  consol*  that  is  connected  to  a  Hil-Std-1750  coeputer 
PASSED  VARIABLES: 

LENGTH-3  -  this  uariabl*  contains  a  count  of  tha  nueber 


RETURNS: 


charge tars  to  print 

NESSHGE-3  -  this  Is  a  location  polntar  for  tha  string  to 
printad 


prints  aassagas  on  user  consol* 

NODULES  CALLED: 

AUTHOR:  Copt.  Steven  A.  Hotchkiss  and 

Capt.  Janrti  far  Fried 

HISTORV:  This  project  was  undertaken  as  a  thesis  project  for 

partial  fulfil  leant  of  requireeents  for  an  NS  degree 
in  Infareation  Science  fro*  tha  Air  Force  Institute 
of  Technology.  Sponsoring  organization  is  the  ASD 
Language  Control  Branch,  Uright  Patterson  RFB.Oh. 


be 


$  4-SEP-85/ 16:09.29  $ 
PRINTOFF 


DO  NOT  LIST  HETAS 


START  OF  META  DEFINITIONS 


DATR8 

LF<0) 

NORNR 
TEST 
_ DAT  AS 


* 

SECTION 


NETA 

3  .  REPEATED  PRESET  NETA 

EQU 

$ 

LOOP 

2, 1,NUN<GF)-I 

void 

GF<_,  1),  NORNR,  _DRTRS 

GOTO 

LABEL 

TEST 

ORTA 

LOOPTEST 

GF<_) 

NETA 

3 

LOOP 

t,  1,GF<_,  1> 

DATA 

LOOPTEST 

MEND 

0F<_) 

LENGTH 

25,9999 

NETA 

0  .  CSECT  NETA 

LOOP 

CSECT 

2,1,31 

SC<0L  **•<_> 


•  GENERATE  REG  EQUATES  META 

* 

REG  NETR 

XNC  LOOP  0, 1, 13 

IC<  R>  EQU  XNC 

LOOPTEST 
NEND 

« 

•  END  OF  NETR  DEFINITIONS 


BASE  FED  EQUATES 


B12 

EQU 

12 

B13 

EQU 

13 

B14 

EQU 

14 

BIS 

EQU 

15 

* 

*  CONDITION  CODE  EQUATES 

* 

-NOP 

EQU 

0 

_LT 

EQU 

1 

-EQ 

EQU 

2 

_LE 

EQU 

3 

jGT 

EQU 

4 

JC 

EQU 

3 

jGE 

EQU 

0 

JCV 

EQU 

8 

-CUT 

EQU 

9 

JCEQ 

EQU 

10 

JCLE 

EQU 

11 

-COT 

EQU 

12 

JCHE 

EQU 

13 

JCGE 

EQU 

14 

JUN 

EQU 

13 

* 

♦  END  OF  EQUATES 

* 

REG 

SECTION 

PRINT 

DEFINE 

PRINTC 

PSSOATRS 

EQU 

3 

PSSCONSt 

EQU 

4 

PSSCOOEI 

EQU 

2 

*  NO  REF  DATA  DECLARATIONS 

*  NO  BVREF/TVPE /ABSOLUTE  DECLARATIONS 

*  LOCAL  AUTONATIC  ORTA  ***  SIZE  IN  UOAOS  —  2  DECIMAL 

*  LOCAL  RUTONRTIC  DATA  FOR  PfiOC  PRINTC 

*  STACK  FRAME  ***  SIZE  IN  UOAOS  —  2  DECIMAL  :  2  HEX 

BK—D03EF  EQU  HEX<0)  .  SIZE 

LENGTH-3  EQU  HEX<0>  .  SIZE 

MESSAGE-3  EQU  HEX<1>  .  SIZE 

*  EN)  OF  LOCAL  AUTONATIC  DECLARATIONS 

*  PSECT  $OATA  IS  EMPTV 


I 


PSSCODE 

PRINTC 


* 

OUTPUT 


* 

LB-J0002 


R2  -  NUMBER  OF  CHRRRCTER8  IN  STRING 
R3  >  LOCATION  OF  CHARACTER  STRING 


ORIGIN  HEX<0> 

EQU  $ 

ORIGIN  HEX (0002 ) 


AISP 

SAR 

BLE 

EQU 

XIO 

TBR 

BEZ 

L 

XIO 

AISP 

SOJ 


R2, 1  .  ADJUST  CHARACTER  COUNT 

R2, 1  .  1ST  TUO  COMMANDS  EQUIVALENT  TO 

R2  ■  A0UNDCR2/2) 

LB— 0002  .  BRANCH  OUT  IF  ILLEGAL  CHAR  COUNT 

$ 

R3,RCS  .  READ  CONSOLE  STATUS 

1,R5  .  CHECK  STATUS  BIT  1 

OUTPUT  .  IP  OFF,  LOOP  BACK  UNTIL  CONSOLE  I 

R5,0,R3  .  GET  NEXT  TUO  CHARACTERS  OF  MESSAGE 

R3,C0  .  PRINT  BOTH  CHARACTERS 

R3, 1  .  POINT  TO  NEXT  TUO  CHARACTERS 

R2, OUTPUT  .  DECREMENT  LOOP  COUNT,  GO  BACK  IF 


EQU  $ 

AISP  R19,2 

POPM  R2,R3 

URS  R1S 

ORIGIN  HEX(OOOO) 

PSHM  R2,R3 

SISP  A 15, 2 

ORIGIN  HEX(0013> 

END 


i 

I 
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tO  October  ISOS 

1.0 

Refhat 

1 

This  routine  is  used  to  convert  ITS  LINK  files  into 
a  format  that  can  be  loaded  into  the  SPERRV  1631 
coeputer  < 1750R  architecture).  The  ITS  files  are 
1  SO'  f i les  and  oust  be  in  the  80  column  record  format 
described  in  the  ENRD  ITS  Load  Nodule  ICO  (COOL  *1009 
contract  •F33657-83-C-0244).  Use  of  the  command  file 
LINK 1750. CON  to  link  all  compiled  modules  mill  insure 
that  these  records  are  of  the  right  format.  The  format 
of  the  SPERRV  loader  records  are  defined  In  Appendix  B 
of  its  programmer  reference  manual.  The  bytes  of  all 
binary  data  fields  must  be  smapped  (i.e.  the  hl<fi 
order  bits  of  a  word  are  smapped  mith  the  lorn  order  8) 
The  only  type  ITS  records  converted  are  binary  and 
end  record  types,  it  also  ignores  all  protection 
indicators,  and  can  not  handle  expanded  memory  jobs. 
Uhen  all  abject  files  are  copied  into  a  single  object 
for  linking  by  the  ITS  LINKER,  the  main  procedure  must 
be  copied  into  the  file  firstlllllll  Otharmise,  this 
application  mill  have  no  may  of  determining  the  point 
that  execution  is  to  begin.  The  'end*  record  created 
by  the  ITS  I  inker  contains  the  lomest  address  of  the 
load  module,  and  this  application  assumes  that  the 
routine  begins  at  that  point.  The  ITS  file  contains 
datafields  that  are  in  HEX  character  representation, 
and  the  SPERRV  1631  expects  binary  data  fields; 
therefore  the  ITS  data  must  also  be  converted  to 
binary 


PASSED  UARIABLES: 

N/fl 

RETURNS: 

N/R 

NODULES  CALLED: 

GetHdr 

Readf 

Printf 

ClrtUp 

IntFil 

MriteRcd 

AUTHOR:  Copt.  Steven  R.  Hotchkiss  and 

_  Copt,  Jennifer  Fried 

H I STORY :  This  project  mas  undertaken  as  a  thesis  project  for 

partial  fulfillment  of  requirements  for  an  NS  degree 
In  Information  Science  from  the  Rir  Force  Institute 
of  Technology.  Sponsoring  organization  is  the  ASD 
Language  Control  Branch,  Uright  Patterson  AFB,0h. 


DATE: 

VERSION: 

rtm. 

NODULE  NUMBER: 
DESCRIPTION: 


START 

ICQNPOOL  < ' loData* ); 
IC0NP00L  (' loCai Is* ); 
ICOMPOOL  < 'RfHtCpI ’ ); 
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PROGRRN  RafNat; 


BEGIN 

ChkSua  •  48*0000'; 

FirstPass  »  Trua; 

LdPt  -  -1; 

Eof  ■  Falsa; 

Buff  »  0; 

BufPtr<0  >  ■  1; 

BufPtr< 1 >  -  1; 

"  Initialize  10  Filas  “ 

IntFI I; 

”  Gat  Haadar  Info  for  Loodar  FI  la  * 

GatHdr; 

UHILE  NOT  Eof; 

BEGIN 

"  Raod  tha  first  80  coluan  racord  ’ 

Raadf( : I tsRcd.Eof ); 

*  Put  Loodar  Info  into  Contiguous  Haaory  Locations  ‘ 
UdsInRcd  ■  CntKO)  -  AsciiO; 

RddrCCO)  -  fiddr<0>; 

Ud1C<0>  -  UdlCO); 

Ud2C<0)  -  Ud2<0); 

UdX(0>  ■  Ud3<0>; 

Ud4C<0>  «  Ud4<0); 

UdSC(O)  *  Ud5<0>; 

UdCC<0)  ■  UdG(O); 

Ud7C(0)  ■  Ud?<0>; 

"  Initial iza  tha  Output  Buffars  * 

FOR  lx:  1  BV  1  UHILE  lx<33; 

CharToBinOx)  -  0; 

FOR  lx:  0  BV  1  UHILE  lx<63; 

OutBuff(lx)  ■  0; 


*  Convert  Char  To  Bin  and  Pack  it  “ 

FOR  lx:  1  BV  1  UHILE  lx<-32; 

BEGIN 

IF  <Rsci  iO<-CharToBinOx))  AND  <CharToBin<  IxX-flsci  i9); 

CharToBinOx)  -  ChcrToBinOx)  -  AsciiO; 

ELSE 

IF  (Rsci  IA<«CharToBlnOx))  AND  <CharToBlnOxX-flsei  IF); 
CharToBinOx)  •  CharToBinOx)  -  AsciiA  +  10; 
HalfBytaOx)  -  NibblasOx); 

END 


IF  Typ<0>  «  *  *; 

BEGIN  "This  Is  a  binary  racord’ 

IF  BufPtr<Buff >  +  UdsInRcd  <«C1  AND  LdPT  -  Laddr<0>; 

BEGIN  "  Old  Racord  and  still  rooa  for  aora  data  fialds 


*  Flip  Flop  the  position  of  eoch  byte  of  a  1750ft  word  * 

FOR  lx:  0  BV  1  UHILE  lx<UdslnRcd; 

BEGIN 

BufByteO(BufPtr<Buff >+lx)  -  FieldLCIx+l); 
BufByteKBufPtr<Buff>+lx>  -  FieldH<lx+1>; 

END 

"  Point  to  eher*  info  fro*  next  ITS  80  coluan  record  “ 

*  is  to  b*  placed  into  this  loadar  r*eord 
BufPtr<Buff>  ■  BufPtr(Buff)  +  Udslnftcd; 

"  Update  load  point  so  th*  next  ITS  racord  can  b*  checked  to  * 
“  sat  If  it  belongs  in  this  loader  record 
LdPt  -  LdPt  ♦  UdsInRcd; 

IF  BufPir(Suff)  -61; 

BEGIN  "  Loader  Record  is  full  and  needs  to  be  aritten  " 

UdsInBuffer  ■  60; 

RcdTypl  ■  ftsci i B; 

UrlteRcd; 

END 

END 


BEGIN  *  Old  record  and  not  enough  rooe  —  or  nee  record  " 

IF  LdPt  -  Loddr<0); 

BEGIN  "  Sane  loader  record,  but  not  enough  room  for  all  ’ 
"  data  fields  in  ITS  record 

"  Seap  Bytes  of  eords  going  into  loader  record  “ 

FOR  lx:  0  BV  1  UHILE  BufPtr<Buff Hlx  <  61; 

BEGIN 

BufByteQ(BufPtr<Buff)+lx>  ■  FieldL<lx+1>; 
BufByte1<BufPtr<Buff>+lx>  -  FieldH<lx+1>; 

LdPt  -  LdPt  ♦  1; 


"  erite  the  full  record  out  " 

UdsInBuffer  ■  60; 

RcdTypI  *  ftsci IB; 

UrlteRcd; 

“  Set  the  load  point  for  this  nee  loader  record  * 

Ldftd(O)  -  LdPt; 

"  Seap  bytes  of  the  other  ITS  data  fields  and  place  the*  into 

*  record.  If  the  next  ITS  record  doesn't  have  the  load  point 

*  coeputed  here,  it  should  be  the  first  entries  for  another 
"  loader  record 

FOR  ly:  lx  BV  1  UHILE  ly  <  UdsInRcd; 

BEGIN 

BufByteO<BufPtr<BuffHly>  *  FleldLOy +  1>; 
BufByte1<BufPtr<Buff Hly)  ■  FieldH< ly+1 ); 

LdPt  -  LdPt  +  1; 


y~ymywFwywyrxwywm)rms'*?rm-'rmy*ymy*J*y'y*?*Jl  v  *.<  mj>?t  -.v  -.» ~.*m'r?\rrTr*m:-r.- 


SO  "Scm  record  not  onogh  rooo" 

ELSE 

BEGIN  11  this  is  ths  start  of  a  naa  loader  record  " 

IF  NOT  FirstPass; 

BEGIN 

IF  BufPtr(Buff)  <>  1; 

BEGIN  "ths  last  record  didn't  gat  fillad  up,  so  it 
"  hasn't  boon  arittan  gat.  Tha  routina 
*  UritaRcd  sats  BufPtr  to  1  bafora  axlt  * 

RcdTgpl  ■  RsciiB; 

UdsInBuffar  *  BufPtr (Buff)  -1; 

UritaRcd; 

END 

END  “  and  not  first  pass  “ 

FirstPass  ■  Falsa; 

“  Sat  tha  load  point  for  this  loadar  racord  * 

Ldfld(O)  -  Laddr(O); 

•  Soap  bgtas  of  IT8  data  fialds  going  into  loadar  racord 
FOR  lx:  0  BV  1  WILE  lx  <  UdsInRcd; 

BEGIN 

BufBgtaO(lx+1)  ■  FialdL< lx+1>; 

BufBgtaKlx+1)  -  FialdHC lx+1 >; 

END 

BufPtr(Buff)  -  UdsInRcd  +  1; 

LdPt  -  Laddr(O)  +  UdsInRcd; 

END  "  and  naa  racord  ” 

END  "and  of  old  racord  not  anough  rooa  —  or  naa  racord  " 
END  ‘and  of  this  is  a  binary  racord" 

ELSE 

BEGIN  ‘this  is  an  axacution  addrass  racord  ” 

IF  Tgp<0)  ■  'E'; 

BEGIN 

RcdTgpl  *  8261;  "blar*  E" 

OutBuff(O)  -  Laddr(O); 

OutBuffd)  *  30;  "  ascii  racord  saparator  " 

UritaRcd; 

END 

END  "and  axacution  addrass  racord  ‘ 

END  ‘and  ahila  loop  " 

"  Urita  and  of  fila  loadar  racord  ” 


RcdTgpl  »  8262;  "  blank  F 


DATE:  10  October  1985 

UERSION:  1.0 

NAME:  Wri tefled 

nOOULE  NUMBER:  7 

DESCRIPTION: 

This  routine  Is  called  by  Refllat  to  do  10  stuff  that 
needs  to  be  done  throughout  the  eain  procedure.  Three 
types  of  SPERRY  1631  loader  records  are  written: 

Binary  , Execution,  aid  End  of  file.  If  the  record  type 
is  a  binary  record,  this  routine  coeputes  a  checksue 
for  it  and  teks  it  on  to  the  end  of  the  record.  Then 
the  record  type  is  art t ten  out  followed  by  the  binary 
record.  If  the  record  type  is  an  execution  record  or 
an  end  of  file  record,  the  record  type  Is  written  out 
followed  by  the  record.  The  variable  ‘Buff  is  a 
global  variable  that  points  to  the  record  to  be 
wri tten. 

PRSSED  VARIABLES:  None 


RETURNS: 
NOOULES  CflLLED: 
AUTHOR: 

HISTORY: 


Nothing 

Printf  —  a  FORTRAN  10  routine 
Capt.  Steven  R.  Hotchkiss  and 
Copt.  Jennifer  Fried 

This  project  was  undertaken  as  a  thesis  project  for 
partial  fulfil  leant  of  requirements  for  an  MS  degree 
in  Inforeation  Science  froe  the  Air  Force  Institute 
of  Technology.  Sponsoring  organization  is  the  ASO 
Language  Control  Branch,  Wright  Patterson  AFB,0h. 


B***» 


I  C0MP00L<  ’  RfMtCp  I  *  ); 

IC0MP00L  ( ’ loOata' >j 

REF  PAX  Printf<RcdTyp, Buffer); 
(LINKAGE  FORTRAN; 

BEGIN 

ITEM  RcdTyp  S  13; 

ITEM  Buffer  C  126; 

END 

DEF  PAOC  UriteRcd; 

BEGIN 

LoopCnt  ■  WdsinBuffer; 

If  AcdTypi  ■  Asci iB; 

BEGIN 

ChkSue  -  48*0000'; 

0utBuff<0)  -  LdAd(O); 

ChkSue  -  ChkSue  XOR  OutBuffB(O); 


OutBuff< 1 )  ■  WdsinBuffer; 

ChkSue  ■  ChkSue  XOR  OutBuf fB< 1 ) 


FOR  lx:  1  BV  1  WHILE  lx  <-  LoopCnt; 
BEGIN 

OUTBUFF < I x+ 1 )  -  BufUd(lx); 

ChkSua  >  ChkSua  XOR  OutBuffB< lx+1 > 
END 

OutBuffB< lx+1 >  ■  ChkSua; 

Pr i ntf <RcdTyp I , OutF I d >; 

END 

ELSE 

Prlntf<RcdTypl  ,0utFld>; 


BufPtr<Buff>  -  1; 
Buff  *  flBS< 1-Buff >: 


DOTE: 

VERSION: 

NRIC: 

NOOULE  NUMBER: 
DESCRIPTION: 


10  October  1989 

1.0 

ioCal Is 
9 


This  coapool  is  required  for  RefNat  to  reference  Its 
associated  FORTRRN  10  routines 


PRSSED  URRIRBLES: 
RETURNS: 

NODULES  CflLLED: 
RUTHOR: 

HISTORV: 


N/fl 

N/fl 

N/R 

Capt.  Steven  R.  Hotchkiss  and 
Capt.  Jenni fer  Fried 

This  project  eas  undertaken  as  a  thesis  project  for 
partial  fulfil  leant  of  requlreeents  for  an  NS  degree 
in  Inforeation  Science  froe  the  Rir  Force  Institute 


of  Technology.  Sponsoring  organization  is  the  RSO 
Language  Control  Branch,  Wright  Patterson  RFB,0h. 


START 

C0MP00L  loCal Is; 

REF  PROC  Wri teficd; 

BEGIN 

END 

REF  PROC  OetHdr; 

ILINKRGE  FORTRRN; 

BEGIN 

END 

REF  PROC  Readft : I tsRcd,Eof ); 
ILINKRGE  FORTRRN; 

BEGIN 

ITEM  ItsRcd  C  80; 

ITEH  Eof  B  1; 

END 

REF  PROC  Printf<RcdTyp,  Buffer); 
ILINKRGE  FORTRRN; 

BEGIN 

ITEM  RcdTyp  S  15; 

I TEN  Buffer  C  126; 

END 

REF  PROC  ClhUp; 

ILINKRGE  FORTRRN; 

BEGIN 

END 

REF  PROC  IntFil; 

ILINKRGE  FORTRRN; 

BEGIN 

END 

TERN 
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DATE: 

UERSION: 

ivifc: 

MODULE  NUMBER: 
DESCRIPTION: 


10  October  1985 

1.0 

loOota 

8 


This  coepool  defines  all  data  required  for  the  JOUIRL 
routine  RefMat  and  its  associated  FORTRAN  10  routines 


PASSED  UARIRBUES: 
RETURNS: 

MODULES  CALLED: 
AUTHOR: 

HISTORV: 


N/A 

M/A 

N/A 

Capt.  Steven  A.  Hotchkiss  and 
Capt.  Jennifer  Fried 

This  project  was  undertaken  as  a  thesis  project  for 
partial  fulfil  leant  of  requlraeents  for  an  MS  degree 
in  Information  Science  froe  the  Air  Force  institute 


of  Technology.  Sponsoring  organization  is  the  ASD 
Language  Control  Branch,  Wright  Patterson  AFB,0h. 


START 


COMPOOL  foOata; 

OEF  ITEM  Infil  C  10 

OEF  ITEM  Outfit  C  10 

OEF  ITEM  FI  I  nan  C  6 

OEF  ITEM  Header  C  80 


OEF  TABLE  ItsTable(O)  U  20; 


BEGIN 

ITEM 

Addr 

C 

4 

POS< 16,00) 

ITEM 

Typ 

C 

1 

P0S< 16,01) 

ITEM 

Cnt 

C 

1 

POS<24,01) 

ITEM 

Cntl 

S 

7 

POS<24,01> 

ITEM 

Udl 

C 

4 

POS<08,03) 

ITEM 

Ud2 

C 

4 

P0S< 16,05) 

ITEM 

Ud3 

C 

4 

P0S<24,0?) 

ITEM 

Ud4 

C 

4 

P0S<00, 10) 

ITEM 

Ud5 

C 

4 

P0S<08, 12) 

ITEM 

Ud6 

C 

4 

P0S< 18,14) 

ITEM 

Ud7 

C 

4 

P0S<24, 16) 

END 

ITEM  ItsAcd 

C 

80; 

OUERLRV  ItsAcd:  ItsTable; 


OEF  TABLE  OutAcd  <0:62)  T  16  U; 
BEGIN 

ITEM  OutBuff  S  15  POSCO, 0); 

ITEM  OutBuffB  B  16  POSCO, 0); 
END 

OEF  ITEM  OutFId  C  128; 

OUERLAV  OutAcd:  OutFId; 


OEF  ITEM  Eof  8  1; 

OEF  ITEM  RcdTypI  S  15; 


OEF  ITEM  RcdTyp  C  2; 

OUERLRV  RcdTyp I :  RcdTyp; 

OUERLRV  Inf i I ,  OutFil,  Filnoa,  Header,  ItsRed,  OutFId,  Eof,  RcdTyp; 


TERR 

i 


I 


I  %• 

i 
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i 


DOTE: 

UERSIQN: 

ill  L. 

riOOULE  NUMBER: 
DESCRIPTION: 


10  Oe lobar  1089 
1.0 

RfMtCpI 

10 

This  cospool  contains  all  ths  variablas  and  tab I as 


PASSED  VARIABLES 
RETURNS: 

MOOULES  CALLED: 
AUTHOR: 

HISTORV: 


that  «*a  ussd  to  unpack  ITS  i inksr  racords,  packs  than 
and  converts  tha  HEX  charoctars  to  binary  data  fields, 
and  than  pi ocas  than  into  a  SPERRY  1631  loadar  racord 
foreat 
:  N/A 
N/A 
N/A 

Copt.  Stavan  R.  Hotchkiss  and 
Copt.  Jami  far  Friad 

This  pro j act  was  undartakan  as  a  thasis  pro j act  for 
partial  fulfillsant  of  raquirasants  for  an  MS  dag raa 
in  Inforaation  Scianca  froa  tha  Air  Forca  Instituta 
of  Tachnology.  Sponsoring  organization  is  tha  RSD 
Languaga  Control  Branch,  Uright  Pa tt arson  AFB,0h. 


ITEM 

ChkSua  B 

16; 

ITEM 

FirstPass  B 

1; 

ITEM 

LdPt  S 

13; 

ITEM 

Buff  U 

8; 

ITEM 

lx  U 

8; 

ITEM 

ly  U 

8; 

ITEM 

UdsInAcd  S 

19; 

ITEM 

LoopCnt  S 

13; 

ITEM 

UdsInBuffar  S 

13; 

ITEM 

Zaro  STATIC 

C 

1  ■ 

•O' 

ITEM 

Asci iO  STATIC 

S 

7; 

AY  Zaro:  Asci iO; 

ITEM 

Nina  STATIC 

C 

1  ■ 

'9' 

ITEM 

Asci i9  STATIC 

s 

7; 

AY  Nina:  Asci 19; 

ITEM 

AA  STATIC 

C 

1  - 

’A' 

ITEM 

Asci IR  STATIC 

S 

7; 

AY  AA 

Asci iA; 

ITEM 

FF  STATIC 

C 

1  ■ 

•F’ 

ITEM 

Asci iF  STATIC 

S 

7; 

AY  FF 

Asci iF; 

ITEM 

BB  STATIC 

C 

2  ■ 

’  B 

ITEM 

Asci iB  STATIC 

S 

13; 

AY  BB 

Asci iB; 

3 


TABLE  LoodPoint  <0>; 
COIN 

ITEM  Ldfid  S  15; 
iMO 


TABLE  BufStuf  <0:0; 
BEGIN 

ITEM  BufPtr  S  7; 
END 


TABLE  Pockedfcd 
BEOIM 

ITEM  AddrC  C  4 
ITEM  MIC  C  4 
ITEM  M2C  C  4 
ITEM  M3C  C  4 
ITEM  M4C  C  4 
ITEM  M5C  C  4 
ITEM  MGC  C  4 
ITEM  M7C  C  4 
END 


<0>  U  8; 

POSCO,  0) 
POSCO,  1> 
POSCO, 2 > 
P0S<0,3> 
POSCO, 4) 
POSCO, 5) 
POSCO, 6) 
POSCO,?) 


TABLE  CharConuart  Cl: 32  >  T  8  U; 
BEGIN 

ITEM  ChcrToBin  S  7  POSCO, 0); 

ITEM  Nibbles  S  3  P0SC4,0); 

END 

ALAS'  PockadRcd:  CharConvart; 


TABLE  HexBuf  0:32)  T  4  U; 

BEGIN 

ITEM  HalfBgU  8  3  POSCO, 0); 

Elf) 

TABLE  Paklts  <0:7)  T  16  U; 

BEGIN 

ITEM  Laddr  S  15  P0SC0,0); 

END 

TABLE  BinFields  C0:7)  T  16  U; 

BEGIN 

ITEM  Field  S  15  POSCO, 0); 

ITEM  FietdH  S  7  PO8C0,0); 

ITEM  FialdL  S  7  P0S<8,0); 

END 

OUEALAV  HexBuf, lx:  Paklts:  BinFialds; 

T  TABLE  DatFlelds  CO)  U  1; 

BEGIN 

ITEM  BufBytaO  S  7  POSCO, 0); 

ITEM  BufBgta t  S  7  P0S<8,0); 

ITEM  BufM  S  15  P0S<0,0); 


[•] 


oooooooooooooooooooooooooo 


(0  October  1983 

1.0 

IntFil 
6 

This  rout in*  is  cat  tad  by  tha  JOVIAL  rout in*  col  lad 
AafNat.  Its  purposa  is  to  proapt  tha  usar  for  tha 
noaa  of  o  fit*  that  *as  croatad  by  an  ITS  link, 
proapt  tha  usar  for  tha  noaa  of  a  f i la  that  tha 
raforaattad  ITS  fila  is  to  b*  arittan  to,  and  than 
opans  both  fllas.  Tha  input  fila  oust  b*  a  ".SO*  f I  la 
and  tha  output  fila  is  a  ".Of IT"  fila. 

PASSED  VARIABLES:  Nona 
RETURNS:  Nothing 

GLOBAL  VARIABLES:  All  uarlablas  usad  arm  global,  and  how*  baan 
dafinad  in  tha  coaaon  (C0MP00L)  cal  lad  loOata 
NOOULES  CALLED:  Nona 

AUTHOR:  Copt.  Stauan  A.  Hotchkiss  and 

_  Capt.  Jam i  far  Friad 

HISTORV:  This  pro j act  acts  undartakan  as  a  th* sis  projsct  for 

partial  fulfil laant  of  raquiraaants  for  an  NS  dogra* 
in  Infaraation  Scianc*  fro*  tha  Air  Fore*  Institut* 
of  Tachnology.  Sponsoring  organization  is  tha  ASO 
Languoga  Control  Branch,  Uright  Pat tar-son  AFB,0h. 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccocccccccccccccc 


DATE: 

VERSION: 

|JAMC  , 

IVIL. 

NOOULE  NUMBER: 
DESCRIPTION: 


Subroutine  IntFil 

IMPLICIT  INTEGER  <A-Z> 

CHARACTER*  Filnaa 
CHARACTER* 10  Infil,  Outfil 

URITE<*,*>'  Enter  Fila  Noaa  <Nax  6  Chcroctars)  ' 

REAO<*, 10)Fi Inoa 
10  F0RMRT(A6) 

I  -  INDEX<Fi  Inaa, ' . '  )  -  1 
IF  (I.LE.O)  THEN 

I  »  INDEX<Fi Inoa, ’  ')  -  1 
IF  (I.LE.O)  THEN 
I  ■  6 
END  IF 
END  IF 

Infil  ■  Fi lnaa(1: I )//' .SO’ 

Outfil  ■  Fi lnaa< 1: I >//’ .DAT' 

WRITE(*,*)‘ Input  Fila  ■  *, Inf  1 1, 'Output  file  ■  *, Outfil 

0PEN<UNIT  -  2,  NAME  -  Infil,  TVPE  •  ’OLD’,  FORM  -  •FORMATTED’ > 
OPEN  (UNIT  ■  3,  NAME  ■  OutFil,  TVPE  ■  'NEW, 

1F0AM  ■  ’ UNFORMATTED •) 


oooooooooooooooooooooooooo 


oooooooooooooooooooooo 


oooooooooonooooooooooooo 


10  October  1985 

1.0 

Printf 

4 

This  routine  is  called  bg  the  JOVIAL  routine  called 
AefNat.  It  Is  used  to  erfte  SPERRV  loader  records  out 
to  a  ” .DAT”  file.  The  name  of  the  fil  being  eritten  is 
stored  in  the  global  variable  Outfit  ehich  eas  set  in 
the  routine  called  IntFil 
PASSED  VARIABLES:  hone 
RETURNS:  Nothing 

GLOBAL  VARIABLES:  All  variables  used  are  global,  and  are  defined  in 
the  coeeon  (COMPOOL)  called  loOata 
hOOULES  CALLED:  None 

AUTHOR:  Capt.  Steven  R.  Hotchkiss  and 

_  Capt.  Jennifer  Fried 

H I  STORY :  This  project  eas  undertaken  as  a  thesis  project  for 

partial  fulfilleent  of  requirements  for  an  US  degree 
in  Information  Science  from  the  Air  Force  Institute 
of  Technology.  Sponsoring  organization  is  the  ASD 
Language  Control  Branch,  Uright  Patterson  AFB,0h. 


DATE: 

VERSION: 
nmt . 

NOOULE  NUMBER: 
DESCRIPTION: 


Subroutine  Printf<RcdTyp,  OutFid) 

IMPLICIT  INTEGER  <R-Z) 
CHARACTERS  RcdTyp 

INTEGERS  0utFld<1:63> 


URITE<3)RcdTyp,<0utFld<l ),  I  ■  1,63) 
HRITE<*,*)'Urite  next  record* 


oooooooooooooooooooooooo 


c 

C  OfiTE: 

C  VERSION: 

C  NRNE: 

C  NODULE  NUMBER: 
C  DESCRIPTION: 


C  PASSED  VARIABLES: 
C  RETURNS: 

C  GLOBAL  VARIABLES: 
C 

C  NOOULES  CALLED: 

C  AUTHOR: 


10  October  IASS 
1.0 
Aeodf 
3 

This  routine  is  called  by  the  JOVIAL  routine  called 
AefMat.  Its  is  used  to  read  80  coluen  records  created 
by  Vie  ITS  linker.  The  none  of  the  fil  being  read  is 
stored  in  the  global  variable  Inf 1 1  ehlch  eas  set  in 
the  routine  called  IntFil.  This  file  Is  a  ".SO"  file 


Nothing 

All  varieties  used  are  global,  and  are  defined  in 
the  coeeon  (C0MP00L)  called  loOata 


C  HISTORV: 
C 


Copt.  Steven  A.  Hotchkiss  and 
Capt.  Jennifer  Fried 

This  project  eas  undertaken  as  a  thesis  project  for 
partial  fulfil  leant  of  requireoants  for  an  NS  degree 
in  Inforeation  Science  from  the  Air  Force  Institute 
of  Technology.  Sponsoring  organization  Is  the  ASO 
Language  Control  Branch,  Wright  Patterson  AFB,0h. 


aooonnooonoooooooonooooo 


ooooooooooooooooooooo 


DOTE:  10  October  1989 

VERSION:  1.0 

NAME:  ClnUp 

MODULE  HUMBER:  5 

DESCRIPTION: 

This  routine  is  called  by  the  routine  RefMat  to  close 
the  files  it  used  for  10 
PASSED  VARIABLES:  None 

RETURNS:  Nothing 

MODULES  CALLED:  None 

GLOBAL  VARIABLES:  Al (variables  used  are  global,  end  are  defined  in 
the  coeeon  (COMPOOL)  called  loOata 
AUTHOR:  Capt.  Steven  R.  Hotchkiss  and 

_  Capt.  Jennifer  Fried 

HISTORV:  This  project  eas  undertaken  as  a  thesis  project  far 

partial  fulfil  leant  of  requirements  for  an  NS  degree 
in  Information  Science  from  the  Air  Force  Institute 
of  Technology.  Sponsoring  organization  is  the  RSO 
Language  Control  Branch,  Uright  Patterson  AFB,0h. 


ooooooooooooooooooooo 


Appendix  D 


96 


I  !  ASNI750  —  Rssaabla  a  (750  sourca  aodula 

;  ! 

i  I  fASNITSO  file 
!  I 

II  fila  ■  input  sourca  naaa  of  aodula  fiia.SI 

I  ! 

I  !  Craata  Assaablar  input  fila  Ul  that  daslgnatas  Mi l-Std-1750fl  as  tha  targat 
I  !  rathar  than  tha  altamata  1750A  targat 

I  ! 

I  CREATE  'Pr.UI 
ASSEMBLE  TARGET-H1750A 
I  TV  'Pr.UI 
I  ! 

I  ASSIGN  'PC. Ul  Ul  !  UPDATE  INPUT  COMMANDS  FILE  (INPUT) 

I  ASSIGN  ‘PC  .SI  SI  I  175QA  ASSEMBLY  SOURCE  FILE  (INPUT) 

t  ASSIGN  ‘PC. OBJ  00  1  OBJECT  OUTPUT 

I  ASSIGN  'PC  .SO  SO  I  SYMBOLIC  OUTPUT 

I  ASSIGN  'PC ,L0  LO  I  LISTING  OUTPUT 

I  ASSIGN  LIB_JOUIfll _ 1730R  01  I  LIBRARY  INPUT 

i  I 

I  SET  VERIFY 
I  M1750A 
i  I 

I  DEASSIGN  SI 
I  DEASSIGN  Ul 
I  DEASSIGN  00 
I  DEASSIGN  SO 
I  DEASSIGN  LO 
I  DEASSIGN  01 
I  ! 

t  DELETE  ’PC  .Ul;* 
t  SET  MOVER  I FY 


J0U1750  —  JOUIflL  COMPILE  FOR  MIL-STO-ITSOR  TRROET 


!  J0U1750  file  I.filetypel  (options] 

! 

!  e.g.,  9J0U1730  TEST1  .SRC  /SVMTRXJOMLV/STflmSTICS 

!  PJ0U1750  TEST2  /MflCHIftXOOE/CROSS 

j 

!  Not*:  If  the  filetype  Is  JOU,  options  say  be  typed  as  2nd  paraeeter. 

I  If  a  filetype  is  suppl ied,  it  eust  be  preceded  by  a  as  shown . 

! 

I  Resulting  object  eodule  has  type  .OBJ 
SET  UERIFV 

JOUIflL  'PI' 'P2' /TflfiGET«1730fl /NO INFO /CROSS /RSSEfl'PS' 


I  LINK  1750  —  Link  one  or  worm  1750A  target  object  Modules. 

! 

!  #L  INK  1750  file 
! 

!  file  ■  object  file  (containing  on*  or  eore  object  Modules) 

!  create  object  file  bg  first  deleting  ail  .obj  files  for 

!  CQNPOOLs  that  don't  contain  any  DEFs.  Then  use  the 

!  following  coee ands  to  create  the  object  file 

I 

!  COPV  *.06J  file.O 
!  RENANE  file.O  file. OBJ 
I 

!  OBJ  files  created  by  the  coepi ler  and  the  asseebler  can  be  copied  to  the 
I  saee  OBJ  file,  but  the  UAX  will  give  an  incoapatible  files  earning.  Ignore 
!  the  earning,  the  copy  is  soda  anyeay 

I 

I  Create  Linker  input  file  “UI" 

j 

SET  UERIFV 
I 

CREATE  'PIVUI 
LINK  DATA,  LI  ST,  DEBUG,  INPUTS 
ALLOCATE  LOCATION- 1000  NODULES  . 

INKEND 

! 


ASSIGN  ‘PI'.UI 

Ul 

!  LINKER  CONTROL  (INPUT) 

ASSIGN  'PI ' .OBJ 

00 

1  OBJECT  MOOULE(S)  (INPUT) 

ASSIGN  'PI' .SO 

SO 

!  LOAD  NOOULE  (OUTPUT) 

ASSIGN  'PI'.LO 

LO 

!  LINKER  LIST  FILE  (OUTPUT) 

ASSIGN  LIB-JOUIAI _ 1750A 

01 

!  LIBRARV  OBJECT  FILE  (INPUT) 

ITSLINK  !  BUN  1750A  Linker. . . reads  logic  device  Ul 

I  Output  on  SO  and  L 0 

DEASSIGN  Ul 
DEASSIGN  00 
DEASSIGN  SO 
DEASSIGN  LO 
DEASSIGN  01 

i 

DELETE  Pr.Ulj* 

SET  NOUEBIFV 


! 

! 


login. con 


This  co— and  procedure  is  invoked  with  each  login, 
and  Mg  be  changed  to  tailor  your  environ— nt. 


!  Set  standard  aliases.  Note  that  several  UNIX-like  aliases  are  set  up 

! 

SET  NOUERIFS' 

SET  PROTECT  I  ON* < SS'STEM : R,0UNER: RUED, GROUP :RU,U0RLD: RUE) /DEFAULT 
SVNBOLS: 

BQ  SHOW  QUEUE /BATCH 

CO  :«•  SET  DEFAULT 

08  DIRECTORS'  /SIZE 

E  EDIT 

HONE  SET  DEFAULT  DSKlflOOL: [ROOL.HOTCHSR 1 

LO  PLOGOUT.CON 

LS  DIRECTORS’ 

PQ  SHOW  QUEUE  SS»S$PfllNT 

PS  :*■  SHOW  PROCESS  !  Like  UNIX  ps  co— and 

PUD  :  —  SHOW  DEFAULT  !  Like  UNIX  ped  co— and 

R  :■*  RUN 

SO  :  —  SHOW  OEUICES 

SO  SHOW  SYMBOLS  /GLOBAL  /ALL 

ST  :  —  SHOW  TERMINRL 

UHO  :*•  SHOW  USERS  !  Like  UNIX  eho  co— and 

SHQ  :  —  SHOU  QUEUE  SLRN$QUEUE/RLL 

S80  SET  TERNINAL/UIDTH-BO 

S132  SET  TERMINRL/UIOTH-132 

J0V1730  PJ0U1730 

LIMKI730:**  PL  INK  1 730 

SIM  1730  PS IM 1730 

ASM  1730  PRSM1730 

UNLOCK  :  —  fUNPROTECT 

! 

I  End  user  defined  keg  ins. 

I 

!  DEFINE  JOUIAL  LIBRRRV  FOR  flUTONRTIC  SEARCHING  FOR  URX  TARGET 

! 

!  ASSIGN  JOUL I BU : JOUL I BU . OLB  LNK$L I  BARRS' 

! 

!  The  fol losing  defines  the  1730R  support  tools  psevdo-co— ends : 

! 

LINK30R  LINKITS 
RRIDX  STOOLS: RAID 

i 

!  END  LOGIN.COM 


•ith  TEXT-10; 
um  TEXT-iO; 

•  ith  TCHEBVSHEF_PflCKflGE; 
um  TCHEBVSHEF-PflCKflGE; 

proc«dur«  TCHEBVSHEF-ECOMOniZflTIOM  is 


This  procedure  is  the  Min  driver  for  ths  Tehebyshef  econosization 
of  a  polynomial . 


ECOMOfl 1 2ED_POLVMOh I flL :  FLOHTJUECTOR  (0. .HflX_DEOflEE)  :■ 

(0.  .(fO-OEOREE  ■>  0.0); 

— The  is  th«  resulting  economized  coefficients  to  ths  polynomial 
SUM:  FLOflTJUECTOR  <0. . ttflXJDEGfiEE )  :«  <0. .HfiXJJEGfiEE  ■>  0.0); 

—This  value  is  a  tsaportry  work  arM  for  ths  sun  of  ths  eolusns 
—  of  ths  work  Mtrix 

U0RICJ1RTRIX:  FLOflTJIRTRIX  <0. .  MRX-DEGREE .  0.  MfiXJDEGREE)  :> 

— Tssporary  work  area  for  forsing  ths  sconosizsd  coefficients 


v  v  v/  ■ 


procedure  DISPLflVJUECTOR  <PRINT_UECTOR:  in  UECTOR)  is 
— Tht  sols  purposs  of  this  rout  ins  is  to  display  an  Intsgsr  wsetor 

packogs  INT-IO  is  nss  INTEGER-10  < intsgsr); 
uss  I NT-10; 

bsgin  — Display  Usetor. 

for  I  in  0. . DEGREE_OF_POLVNON I  PL  loop 
put  <1  >; 
put  <" 

put  < PRINT JJECTOR  (I >>; 
nsst_i  ins; 
snd  loop; 

snd  OISPLflVJJECTOR: 


procedure  DISPLRV-FLOflTJUECTDR  (PRINTJUECTOB:  in  FLOAT  JJECTOR)  Is 
—The  sol*  purpose  of  this  routine  is  to  display  a  floating  point  vector 

package  INT-IO  is  nee  INTEGER.  10  (integer); 
use  INT_I0; 

package  FLT-IO  is  nee  FLOAT- 10  (float); 
use  FLT-IO; 

begin  — Display  Float  Uector. 

for  I  in  0. .DEGAEEjOF-POLVTOTIflL  loop 
put  (I); 
put  C  “>; 

put  (PRINT-UECTOR  (l»; 
nee_l  In*; 
end  loop; 

end  0 1  SPLRV-FLORT  JJECTOR : 


procedure  DISPLflVJIfiTRIX  (PfllNTJIfiTRIX:  in  MBTRIX)  is 
— The  sols  purpose  of  this  routine  Is  to  display  an  integer  eatrlx 


package  INT-IO  is  nee  INTEGER- 10  (integer); 
use  INT-IO; 

package  FLT-10  is  nee  FL0RT-I0  (float); 
use  FLT-10; 

beg  in  — Display  Matrix. 

far  I  in  0. .OEGREE_OF_POLVNOMIflL  loop 
put  (I); 
put  (*  “>; 

for  J  in  0. .DEQF€E_OF_POLVNONIflL  loop 
put  (PRINT-NRTRIX  (|,J)>; 
put  C  “); 
end  loop; 
nee-l  ine; 
end  loop; 

end  OISPLRVJVrmiX; 


» 


procedure  0 1 SPLflV_FLOflT_HflTR I X  (PR I NT-MOTRI X :  in  FLORT-MRTRIX)  is 
—The  sole  purpose  of  this  routine  is  to  display  a  floating  point  eatrix 

package  INT_IO  is  nee  INTEGER_IO  (integer); 
use  I NT-1 0; 

package  FLT-IO  is  nee  FLORT-IO  (float); 
use  FLT-IO; 

begin  — Display  Float  Matrix. 

far  I  in  0. .DEGREEjOF-POLYNONIRL  loop 
put  (I); 
put  (*  “>; 

for  J  in  0. .DEOREE-OF_POLYNOrilflL  loop 
put  ( PRINT JIRTRIX  (l,J)>; 
put  C  "); 
end  loop; 
nee_l  ine; 
end  loop; 

end  OISPLflV_FLOflT_NfiTRIX; 


gin  — Tchebyshef  Econo* i ration. 


INPUT-COEFFICIENTS; 
put  ("Input  Coefficients”); 
im.1  in*; 

0 1 SPUAV-FLOAT JUECTOR  (COEFFICIENTS); 

COMPUTE_TCHEBVSHEF_POLVNOf1 1  flL ; 
new -I  in*; 

put  (“Tchebyshef  Polgno*ial"); 
n**_l  in*; 

DISPLflV-HfiTBIX  (TCHEBVSHEF-POLVNONiflLS); 

CONPUTE_POUERS_OF_TCHEBVSHEF; 
n**i_l  In*, 

put  ("Powers  of  Tchebyshef"); 
n**_t  in*; 

DISPUW-FLOflTJIRTRIX  (POOSjOF-TCHEBVSHEF); 

— Generate  the  work  eatrix  used  in  the  final  calculations  of  the  economized 
—  polynomial.  Again  the  eatrix  is  lower  triangular, 
for  I  in  0. ,DEGAEEjQF_POLVNOniflL  loop 
f or  J  in  0. . I  loop 

UORK-MATRIX  (l,J)  f  I  oat  (MULTIPLIER  (l»  *  POUERSlOF-TCHEBVSHEF  (l,J) 

*  COEFFICIENTS  (I); 
end  loop; 
end  loop; 

— Accumut  *•  the  sue  of  the  work  eatrix  columns 
for  I  in  0. . OEGREE-OF-POLVNON I AL  loop 
for  J  in  0. .DEGPEEJOF-PQLVNOtllAL  loop 
SUN(J)  :■  SUTI(J)  +  (UOAK_NATAIX( I ,  J>>; 
end  loop; 
end  loop; 

— Perform  the  final  additions  and  multiplications  to  form  the  result, 
for'  I  in  0. . (DEGAEEjOF-POLVNONIAL  -  1)  loop 
lor  J  in  0. . I  loop 

ECONOmZEDJPOLVHONIAL  (J)  :■  ECONOniZED-POLStCniAL  (J)  ♦ 
f  1  oat(TCHEBVSHEF_POLVNOM I ALS  (l,J»  *  SUM  (I); 
end  loop; 
end  loop; 

new_ I  In*; 

put  ("Economized  Polynomial'); 
new- 1  ine; 

D I  SPLAY .FLOAT -UECTOA  (ECONOMIZED-POLVNOniflL); 


end  TCHEBVSHEF-ECONOh 1 2AT I  ON : 


DATE:  Dkm bar  1,  1965 
VERSION:  1.0 


STRING-TO-INT 

INPUT-COEFFICIENTS 

COMPUTE-TCHEBVSHEF  -POLYNON I RL 

CONPUTE-POMERSJJF-TCHEBYSHEF 

DISPLAY-VECTOR 

DISPLAY-MATRIX 

DESCRIPTIONS:  Provided  with  aoch  routine. 

PASSED  VARIABLES:  The  input  to  this  systea  is  the  description  of 
the  polynoaial  to  be  econoeized. 

RETURNS:  The  result  of  processing  is  the  coefficients  of  the 
econoeized  oolunoeial . 

CALLING  NOODLES:  TCHEBYSHEF-ECONOMIZATIQM 

AUTHOR:  Capt  Jennifer  Fried  and 
Capt  Steven  Hotchkiss 

HISTORV:  Original  version,  Dec  1,  1965 


eith  TEXT-10;  use  TEXT-10; 
package  TCHEBYSHEF-PACKAGE  is 


—  This  package  receives  the  coefficients  of  a  polynomial  that  is  to  be 

—  econoeized,  computes  its  Tchebyshmg  polynoaial,  and  the  poeers  of 

—  Tchebyshef  eatrix. 


— Unconstrained  type  declarations 

type  MATRIX  is  array  (integer  range  <>,  integer  range  <>)  of  integer; 

— Matrix  of  integer  values,  used  to  contain  the  Tchebyshef  polynomials 
type  FLOAT-MATRIX  is  array  (integer  range  o,  integer  range  <>)  of  float; 


— Matrix  of  floating  point  values,  used  to  contain  the  poeers  of 
—  Tchebyshef 

type  VECTOR  is  array  (integer  range  <>)  of  integer; 

— Vector  of  integer  values,  used  to  contain  the  multiplier  of  the  eatrix 
type  FLOAT-VECTOR  is  array  (integer  range  <>)  of  float; 

— Vector  of  floating  point  values,  used  to  contain  the  coefficients  of 

—  the  polynomial 

— Variable  declarations 
MAX-DIOIT :  integer  :■  19; 

— The  eaxieue  nueber  of  digits  pereitted  in  a  number  is  nine. 

—  This  value  represents  the  eaxieue  input  string  length  for  teo  numbers 

—  and  a  slash,  “/". 

MAX-DEGREE:  integer  :■  9; 

— The  eaxieue  value  of  the  largest  exponent  of  the  polynoaial 
DEGREE-QF-POLVNOMIAL:  integv  :«  0; 

—The  actual  value  of  the  largest  exponent  as  input  by  the  user 
COEFFICIENTS:  FLOAT-VECTOR  (0. .MAX-DEGREE)  :■  (0. .MAX-DEGREE  ->  0.0); 
—Contains  a  coefficient  for  each  degree  of  the  polynoaial  that  was 

—  specified  by  the  user 

MULTIPLIER:  VECTOR  (0. .MAX-DEGREE)  :■  (0. .MAX-DEGREE  ->  0); 

—This  vector  contains  the  reciprocal  of  the  values  contained  on 

—  the  diagonal  of  the  Tchebyshef  polynomial  eatrix. 

—  Used  in  generating  the  econoeized  polynomial. 


TCHEBVSHEF_POLVMOHIflLS:  MfiTRIX  <0. .HRX-DEOREE,  0. .MRX-DEGREE) 

<0.  .IlflXJDEOREE  ->  <0.  .HfiX_DEGflEE  ->  0»; 

—Tha  Mtrix  obtained  ahan  using  tha  Tchabyshaf  foraula. 
POWERSJOF_TCHEBVSHEF :  FLOfiTJIRTRIX  <0. .MRX-DEGREE,  0. . MRX-DEGREE > 

<0. .MRX-DEOREE  »  (0. .MRX-DEGREE  ->  0.0)); 

— ' Tha  Mtrix  foraad  shan  applying  tha  sacond  stop  of  tha  aconoaization 

—  algor i tha 

function  STR1NG-T0-INT  (S:  string)  ratum  intagar; 

— This  function  is  usad  to  convart  tha  input  coafficiant  string  into  cat 

—  intagar  uaiua  that  aquatas  to  tha  nuaarator  and  tha  danoainator. 

— Thasa  procaduras  parfora  tha  functions  spacifiad  by  this  pockoga 
procadura  INPUTjCOEFFICIENTS; 

— Oat  tha  Input  coaff lei ants  for  tha  polynoalal 
procadura  COhPUTE-TCHEBYSREF-POLVNOri  I  PL; 

— Oanarota  tha  Tchabyshaf  polynoalal  Mtrix 
procadura  COIVVrE^OUERSJOP-TCHEBVSHEF; 

—Oanarota  tha  poaars  of  Tchabyshaf  Mtrix 


and  TCHEBVSHEF-PRCKRGE; 


function  STRING-TOUNT  <8:  string)  return  integer  is 
— String  to  integer  equivalent  conversion. 


CHRR 

DIGIT 

MULTIPLIER 

FINRL-RESULT 

POSITION 


:  character;  — Individual  nueber  in  each 

—  placeholder  of  the  input  string. 

:  integer;  —Individual  nueber  in  each  placeholder 
—  of  the  output  integer. 

:  integer  :■  t;  —Tens  value  of  the  integer 
—  pointer. 

:  integer  :■  0;  —Output  integer  being  generated. 

:  integer  :■  S' last;  — Pointer  into  input  string 

—  (eoves  right  to  left). 


Mg in  — String  to  integer  conversion. 

— Starting  froe  the  end  of  the  input  string,  process  each 
—  successive  character  until  all  characters  have  been  converted, 
ehi le  POSITION  »  S'first  loop 

— Get  one  character  digit  froe  the  input  string. 

CHRR  :■  S(POSITION); 

—If  this  is  a  valid  character  digit  representation,  convert  the 

—  character  into  its  nueeric  representation,  and  eultiply  it  by 

—  its  tens  value. 

if  CHRR  in  'O'. .'9'  then 

DIGIT  :■  character'posiCHAR)  -  character 'pos< ‘O' >; 

DIGIT  DIGIT  *  MULTIPLIER; 

—  If  the  final  value  si  1 1  be  the  eost  negative  nueber, 

—  designate  it  as  the  eost  negative  nueber  and  stop 

—  processing.  The  reason  this  is  done  is  to  adjust  for  the 

—  problee  that  the  absolute  value  of  the  eost  negative  nueber 
—  Is  1  digit  larger  than  the  eost  positive  nueber  and  sill 

—  result  in  an  out-of-bound  condition. 

if  integer’ last  -  (FINRL-RESULT  -  1)  +  DIGIT  then 
FINRL-RESULT  integer' first; 

POSITION  :■  S'first; 
else 

— Others ise,  this  is  not  the  eost  negative  nueber.  Thus, 

—  add  the  current  digit  to  the  rest  of  those  found,  and 

—  increment  the  tens  value  to  the  next  larger  nueber. 
FINRL-RESULT  :•  FINRL-RESULT  *  DIGIT; 

MULTIPLIER  :■  MULTIPLIER  *  10; 
end  if; 

—If  the  original  input  eas  negative,  then  negate  the  results, 
elsif  CHRR  -  '-'  then 
FINRL-RESULT  :■  -FINRL-RESULT; 
end  if; 

— Adjust  the  pointer  into  the  input  string  to  point  to  the  next 


procedure  I NPUT  JCOEFF I C I  ENTS  is 

—This  procedure  obtain*  the  information  about  the  input  polynomial  and 

—  convert*  the  coefficient*  into  floating  point  format 

package  I  NT -.10  is  nee  I  NTEGER_I0<  integer); 
u*e  I NT-10; 

POWERS :  integer  :■  3; 

—Indicate*  ehether  all  powers,  only  the  even,  or  only  the  odd  poeer* 

—  are  present  in  the  input  polynomial.  Originally  set  to  out  of 

—  bound*  condition  to  verify  proper  input. 

STEPS:  integer  :■  2; 

— Increment  value  for  entering  the  coefficient*  of  the  polynomial 
INITIAL:  integer  0; 

— Starting  value  for  the  value  of  the  exponent 
COUNTER:  integer; 

— Loop  counter  through  the  input  string 
NUMERATOR:  integer; 

— Numerator  of  the  coefficient 
DENOMINATOR:  integer; 

— Denominator  of  the  coefficient 
CONVERT -STRING:  string  <1. .MAX-DIGIT); 

— String  representation  of  the  coefficient 
LAST-DIGIT :  integer; 

—Actual  length  of  the  input  string 

begin  —Input  Coefficient*. 

— Obtain  the  value  of  the  largest  exponent  of  the  polynomial. 

—  It  must  be  betmeen  2  and  9. 

•hi  le  OEGREEJOF-POLVNOMIAL  <  2  or  DEGREEJOF-POLVNOM I AL  >  MAX-DEGREE  loop 
put  ("Enter  the  degree  of  polynomial  desired.  (Minimum  is  2>:  ">; 

get  (OEGREEJOF-POLVNOMIAL); 
new-l  ine; 
end  loop; 

—Obtain  an  indicator  for  the  type  of  the  polynomial’s  exponents 
•hi le  POWERS  <  0  or  POWERS  >  2  loop 
put  ("Enter  0  for  coefficients  for  ALL  powers  of  X"); 
new- I ine; 

put  ("Enter  1  for  coefficients  for  000  powers  of  X"); 

put  ("Enter  2  for  coefficients  for  EVEN  powers  of  X" ); 
new. I  ine; 
get  (POWERS); 
end  loop; 

—Set  the  initial  and  incremental  values  for  obtaining  the  polynomial 

—  coefficients.  Saves  time, 
if  POWERS  -  0  then 

STEPS  :■  1; 
elsif  POWERS  -  1  then 
INITIAL  :■  1; 
end  If; 


— Obtain  tha  coefficients  for  each  eieeent  of  the  polynomial 
put  ("Enter  the  coefficients  of  the  series  being  expanded  by"); 

put  <“  entering  a  fraction,  i.e.  -2/3  or  +2/3  or  2/3" >; 
nee_l ine; 

put  (“Coefficient  for  X**  ">; 
nee- I ine; 

— Loop  through  all  elements 
ohile  INITIAL  <■  OEGREE-DF-POLVNOM I AL  loop 
put  (INITIAL); 
put  (‘  «  "); 

get-line  (C0NUERT-STRING,LA8T-DI0IT>; 
nee_ I ine; 

COUNTER  1; 

— Step  through  the  Input  string  looking  for  the  '/*  ehlch  separates 

—  the  numerator  froe  the  denominator.  If  one  does  not  exist,  or  ft 

—  appears  in  either  the  first  or  the  last  position  in  the  string, 

—  then  the  coefficient  must  be  reentered, 
while  COUNTER  <«  LAST-DIGIT  loop 

if  (CONUERT-STRINO  (COUNTER)  -  ‘/’)  and 
(COUNTER  /-  CONUERT -STRING' f irst  and 
COUNTER  /»  LAST-DIGIT)  then 

declare 

— Once  the  "/"  has  been  located  and  is  in  a  proper  location 

—  obtain  the  numerator  string  and  the  denoeinator  string. 
NUMERATOR-STRING:  string  renames 

CONUERT-STRINO  (COfk€RT-STRING‘  f  irst. .  (COUNTER  -  1»; 
DENOMINATOR-STRING:  string  renames 
CONUERT-STRINO  ((COUNTER  +  1 ). .LA8TJ3IGIT); 

begin  — Block 

— Convert  the  two  strings  into  integers 
NUMERATOR  STR I NG-TO-I NT  (NUMERRTQR-STRING); 

DENOMINATOR  :■  STR I NG-TO-I  NT  (DENOMINATOR-STRING); 

—  If  the  denominator  is  a  valid  value,  then  generate  the  floating 

—  point  value  for  the  coefficient 
if  DENOMINATOR  /-  0  then 

COEFFICIENTS  (INITIAL)  f loat(NUNERATOR)  /  float (DENOMINATOR) 

— Increment  to  the  next  element  in  the  polynomial. 

INITIAL  :■  INITIAL  +  STEPS; 
end  if; 

—  Indicate  that  this  coefficient  has  been  found  and  converted 
COUNTER  LAST-DIGIT; 
end;  — Block 

end  if; 

— Point  to  the  next  character  in  the  input  string 
COUNTER  :■  COUNTER  +  1; 
end  loop; 
end  loop; 

end  INPUTjCOEFFICIENTS; 


procedure  COOPUTE_TCHEBVSHEF_POLVNOt1 1 flL  is 
— Generate  the  matrix  of  the  Tchebyshef  polynomial .  The  procedure  uses 

—  values  of  the  eatrix  elements  that  have  already  been  found. 

—  The  algorithm  is  recursive  in  *!at  respect. 

begin  — Compute  Tchebyshef  Polynomial . 

— The  first  teo  elements  oust  be  initialized  to  ailoe  the  fol toeing 

—  passes  to  use  them. 

TCHEBVSHEF-POLVNOfl I ALS  (0,0)  1; 

TCHEBVSHEF-P0LVN0HIRL8  (1,1)  :■  1; 

— Loop  through  the  loeer  tricngular  portion  of  the  matrix 

—  and  calculate  the  Tchebyshef  polynomial  values. 

for  I  in  2. .HflX_DEGREE  loop 
for  J  in  0. . I  -2  loop 
TCHEBVSHEF-POLVNOI I flLS  <I,J>  :* 

TCHEBVSHEF-P0LVTK3M I  flLS  <I,J>  -  TCHEBV8HEF -PXVMOH I  flLS  (I  -  2,J) 
end  loop; 

for  J  in  0. . I  -  1  loop 
TCHEBVSHEF_P0LVM0t1 1  flLS  <I,J  +  1)  :« 

TCHEBVSHEF-POLVMOn I flLS  <I,J  +  1)  ♦ 

<2  *  TCHEBVSHEF_P0LVM0n  I  flLS  (I  -  l,J»; 
end  loop; 
end  loop; 

end  COOPUTE_TCHEBVSHEF_POL VNOM I flL : 


procedure  COfFLTE-POUERS-OF.TCHEBVSHEF  is 
— Coepute  the  eatrix  for  the  poeere  of  Tchebgshef 

COEFFICIENT-LIST:  FL0flT_t€CT0R  (0. .MRX-DEGREE)  :» 

<0. . MRX-DEGREE  »  0.0); 

INDEX:  integer  :»  OEGREEJOF-POLYNON I flL ; 

STEP:  integer; 

POINTER:  integer; 

begin  — Coepute  Powers  of  Tchebgshef. 

•hi le  INDEX  >■  0  loop 
MULTIPLIER  (INDEX)  :•  1  /  TCHEBVSHEF-POLVNON I RLS  (INDEX, INDEX); 

STEP  :-  INDEX; 

•hi le  STEP  >*  0  loop 

COEFFICIENT-LIST  (STEP)  :-  f I  oat (TCHEBVSHEF-POLVNON I RLS  (INDEX,STI 
STEP  :-  STEP  -  1; 
end  loop; 

POUERS-OF-TCHEBVSHEF  ( INDEX, INDEX)  :■  f.O; 

STEP  :■  INDEX  -  2; 

•hi le  STEP  >«  0  loop 
POUERS-DF-TCHEBVSHEF  (INDEX, STEP)  :- 
-  (COEFFICIENT-LIST  (STEP)) 

/  f I  oat (TCHEBVSHEF-POLVNON I RLS  (STEP, STEP)); 

POINTER  STEP; 

•hi le  POINTER  >■  0  loop 
COEFFICIENT-LIST  (POINTER)  :■ 

COEFFICIENT-LIST  (POINTER)  ♦  POMERSJOF-TOCBVSHEF  (INDEX, ST 
*  float(TCHEBVSHEF_POLVNOff/flLS  (STEP, POINTER »; 

POINTER  :-  POINTER  -  2; 
end  loop; 

STEP  :■  STEP  -  2; 
end  loop; 

INDEX  :■  INDEX  -  1; 
end  loop; 

end  CONPUTE-POUERSJOF-TCHEBVSHEF , 


end  TCHEBVSHEF.PRCKRGE; 


—  Oate: 

I  —  Version: 

—  Ncmw: 

—  Module  Nueber: 

—  Description: 

i 

—  Passed  Variables: 
—  Returns: 

—  Globals  Used: 

—  Modules  Cal  led: 

1  —  Author: 

—  History: 


28  Noveeber  1089 
1.0 

flpprox_Dr i ver 
1.0 

This  routine  loops  until  a  user  is  done  approx i eating 
ehicheuer  function  he  desires 


Choice 

MENU 

Copt.  Steven  A.  Hotchkiss  and 

Capt.  Jennifer  Fried 

Developed  as  a  thesis  and  ADR  project 


•ith  GL06AL-DRTABRSE;  use  GLOBAL-DATABASE; 

with  APPROXIMATORS;  use  APPROXIMATORS; 

j  eith  TEXT-10;  use  TEXT_I0; 

procedure  APPROX-DRIVER  is 

MUM:  integer  :■  0; 

DEN:  Integer  :■  1; 

i  CHOICE,  KEY  :  character; 

QUIT  :  character  :-  '7'; 

package  IMT_I0  Is  nee  I NTEGER-I0< INTEGER); 
use  I  NT-10; 

package  FLT-IO  is  nee  FLORT_IO<LOHG_FLORT>; 
use  FLT-IO; 

begin 

set-page-1  ength<24  >; 

—  initialize  data  points 
COMPUTE-TCHEBYSHEV; 

—  let  the  user  approxleate  as  eany  functions  as  needed 


6 


■hil«  (CHOICE  /»  QUIT)  loop 


—  smimct  function  to  approximate 

—  by  giving  users  a  Mnu  of  options 
HEMK CHOICE); 

—  ust  the  built  functions  to  make  a  more  accurate  approximation 

rar»vrE_pf«6_flpp^ 

COnPUTE-IfK; 

if  CHOICE  /•  QUIT  than 
for  I  in  0. .H  loop 

if  C(NUf1,l>  /■  0.0  or  C(Dmn, I )  /■  0.0  then 
put (“a**"); 
put<  I); 
putr  *■>  "); 
put(C(HUn, I )>; 
put<“ 
put<  Im¬ 
pute  *»>  “); 
put(C(0EN, I )>; 
nmm_l  inm; 
mnd  If; 
mnd  loop; 
mnd  if; 

putCHit  any  kmy  to  continue"); 

gmt(KEV); 

nmm_l inm; 


mnd  loop; 


mnd  flPPROOCJQRIUER; 


—  Data: 

—  Version: 

—  Nose: 

—  Module  Nueber : 

—  Description: 

—  Passed  Variables: 

—  Returns: 

—  Globais  Used: 

—  Modules  Called: 

—  Author: 

—  History: 


28  Novesber  1985 

1.0 

GL08RL-DRTRBR8E 

2.0 

Contains  all  global  uariables 

M/fl 

M/fl 

ALL 

H/A 

Copt.  Steuen  A.  Hotchkiss  and 

Capt.  Jennifer  Fried 

Coepleted  for  Thesis  and  ADA  project 


package  GLOBAL-DATABASE  is 


type  LONG-FLOAT  is  digits  9; 

type  VECTOR  is  array( integer  range  0..25)  of  LONG-FLOAT; 
type  MATRIX  is  array< integer  range  0..25,  Integer  range  0..25)  of 
LONG-FLOAT; 

type  PAOE-MRTRIX  is  array  (integer  range  0. .25,  integer  range  0. . 1, 

integer  range  0. .25)  of  LONG-FLOAT; 


T:  MATRIX; 

A:  PAOE-MRTRIX; 


D:  VECTOR; 

M:  Integer; 

K:  integer; 

M:  integer; 
MACLAURIN:  VECTOR; 


COEFFICIENT:  string<1.  .33); 

EPS:  LONG-FLOAT; 

C:  MATRIX; 


—  Matrix  containing  the  coefficients  of 

—  different  posers  of  Tchebyshev  polynoeials 

—  Used  to  contain  the  series  of  PACE  approx 

—  R(S,N  or  D,C) 

— ■  S  is  the  series  nueber 

—  M  or  D  N  -  0  for  the  nuserator 

—  0  -  1  for  the  denoeinator 

—  C  -  coefficient  for  a  poser  of  X  for  the 

—  particular  series'  nuserator  or 

—  denoeinator  _ 

—  Error  values  of  PACE  approx I  sat ions 

—  Poser  of  the  nuserator  polynoelal 

—  Poser  of  the  denoeinator  polynoeiai 

—  Poser  of  the  initial  poser  series 

—  Contains  the  coefficients  for  the 

—  different  posers  of  “X"  for  the  poser 

—  series  expansion  of  a  function 

—  Used  to  contain  user  entered  coefficients 

—  a  poser  series  expansion 

—  Convergent  epsl Ion 

—  Final  rational  approxisation 


end  GLOBAL-DATABASE; 


—  DATE: 

—  Uersion: 

—  Mom: 

—  nodule  Nuaber: 

—  Description: 

—  Passed  Variables: 

—  Returns: 

—  Olobals  Used: 

—  Hodules  Cal  led: 

—  Author: 

—  History: 


28  November  1985 

1.0 _ 

connoN-PAOCs 

3.0 

This  package  contains  procedures  that  are  invoked 

throughout  the  systee 

h/A 

N/A 

H/A 

H/A 

Capt.  Steven  A.  Hotchkiss  and 

Copt.  Jennifer  Fried 

Developed  as  a  thesis  and  ADA  project 


with  GLOBAL-DATABASE;  use  0L08ALJDATABASE; 
package  COHMOH-PAOCS  is 

procedure  POUER_PROnPT<HUH,  DEN:  out  integer;  Epsilon:  out  L0HG_FL0AT>; 
procedure  GET_C0EFFICIENT8<STRUCTURE: in  character;  POWER:  in  integer); 
function  PROOUCTtFROH,  TO,  BV:  integer)  return  LOHG-FLOAT; 
function  FACTORIAL  (HUHBER:  integer)  return  LOHG-FLOAT; 
end  COmOH-PROCS; 


1  !9 


■  i-ii.  i  a.  e.  . 


'«  *■ **'  h",  e.'  a  *.  A-'.  iu*.  .eJl 


•ith  TEXT_IO;  un  TEXT-10; 

package  body  COWOH-PROCS  is 


pack cty«  INT-IO  is  now  INTEGER_IO<  integer); 
use  INT-IO; 

package  FLT-IO  is  r*m  FLOflT_IO<LONO_FLOflT>; 
use  FLT-IO; 

proetdurt  POUER-PfK)NPT<NUN,  OEN:  out  integer;  Epsilon:  out  L0NQ_FL0flT>  is 
begin 

set_page_l  ength<24  ); 
loop 

put< "Enter  the  poor  of  the  nu»erator<»ust  bo  integer)  “>; 

got<NUN); 

no«-l ino; 

putCEntor  th#  pooor  of  tho  denoeina tor (oust  bo  intogor)  ">; 

get(OEN); 

noo-l ino; 

put<"Enter  Um  opsilon  of  convergence."); 

putCThis  oust  bo  a  roal  fraction  and  ontorod  as  0.x"); 

put<”Uhora  x  is  any  siring  of  digits  up  to  9  in  length 

got<EP8IL0N>; 

noo_l ino; 

oxit; 
and  loop; 
axcaptian 

ohon  data— arror  *> 

put-1 ine<"lnval id  Entry.  Roontor  data”); 


and  POWER -PR0NPT; 
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procedure  GETjCOEFFICIENTS(STRUCTURE:  in  character;  POWER:  in  integer)  is 


COEFF  :  L0NG_FL0flT  :■  0.0; 

FROM,  TO,  BV  :  integer; 

procedure  GET-POWER (NUMBER:  in  integer;  COEFF:  out  LONG-FLOAT)  is 

LAST -SWAP  :  integer  :•  COEFFICIENT' last-1; 

NUNERflTOR  :  boo I  son  :■  TRUE; 

OUT.COEFF  :  LONG-FLORT; 

CHflR-PTR  :  integer; 

NUN,  DEN  :  stringd. .  15); 

INPUT-ERROR  :  exception; 


procedure  CONPUTE_flERL-COEFF(NUH,  DEN:  in  string; 

COEFF:  out  LONG-FLORT >  is 

CHRR-PTR:  integer  :■  NUN' first; 

NUNERflTOR:  LONG-FLORT  :■  0.0; 

OENONINRTOR,  SIGN:  LONG-FLORT  :-  1.0; 


begin  —  CONPUTE-REfll _ COEFF 

if  (NUH(CHflR_PTR )  -  '  +  ')  then 
CHRR-PTR  :-  CHRR-PTR  ♦  1; 

•Isif  (NUN(CHflR-PTR)  »  then 
SIGN  :■  -SIGN; 

CHRR-PTR  :-  CHRR-PTR  ♦  1; 
end  if; 

while  «NUN(CHflR_PTR)  /■  '  * )  and  <CHflR_PTR  <*  NUN' last))  loop 
NUNERflTOR  :-  NUNERflTOR  *  10.0  + 

LONG-FLOAT (charac ter' pos< NUN <CHflR_PTR)>  - 
character 'pos< 'O' )); 
end  loop; 

CHflR_PTR  :■  DEN' first; 
if  (DEN(CHflR-PTR)  -  '♦' )  then 
CHflR-PTR  :»  CHflR-PTR  +  1; 
elsif  < OEM < CHflR-PTR)  ■  '-')  then 
SIGN  :■  -SIGN; 

CHflR-PTR  :»  CHflR-PTR  +  1; 
end  if; 

while  < (DEN (CHflR-PTR )  /■  •  ' )  and  (CHflR-PTR  <■  DEN' last))  loop 
NUNERflTOR  :-  NUNERflTOR  *  10.0  + 

LONG-FLORT (charac ter 'pos(DEN( CHflR-PTR))  - 
character 'pos( '0' )); 
end  loop; 

COEFF  NUMERATOR /OENON I NflT0R*S I GN; 


end  COfPUTE-PEflL-COEFF ; 


begin  —  GET-POWER 


nMLpcgt) 

loop 

—  proop t  tht  usor 

put ("Enter  the  coefficient*  for  x**“>; 

put(NUMBER); 

putC  — >  ">; 

get (COEFFICIENT); 

—  pock  and  seperate 

for  I  in  COEFFICIENT' range  loop 

if  CO'  <-  COEFFICIENTd)  and  COEFFICIENTd  >  <■  'Q'>  or 
COEFFICIENT(I)  ■  or  COEFFICIENTd)  -  or 
COEFFICIENTd)  -  '/'  or  COEFFICIENTd)  -  '  '  then 

if  COEFFICIENTd)  »  '/*  then 
CHflR-PTR  :■  OEN* first; 

NUMERATOR  :■  FALSE; 

elsif  COEFFICIENTd)  -  '+'  or  COEFFICIENTd)  *  or 

CO'  <-  COEFFICIENTd)  and  COEFFICIENTd)  <-  '9') 

then 

if  NUMERATOR  then 

NUM(CHflR-PTR)  :■  COEFFICIENTd); 

else 

OEN(CHflR-PTR)  :»  COEFFICIENTd); 
end  if; 

CHflR-PTR  :*  CHflR-PTR  ♦  }; 
end  if; 

else 

raise  INPUT-ERROR; 
end  if; 
end  loop; 

exit; 
end  loop; 

COMPUTE-RERL-JCOEFFtNUM,  DEN,  OUT-COEFF); 

COEFF  :«  OUTJCOEFF; 
put (OUTJCOEFF); 
nee_ I ine; 


exception  _ 

•hen  INPUT-ERROR  ■>  put_l  ine<" Input  Error.  Reenter  value."); 
nee_ I ine; 


end  GET-POUER; 


gin  —  GETjCOEFFICIENTS 
se  t_page-l  ength<24 ) ; 


new_poge; 


put_l  ine< "Enter  the  coefficients  for  each"); 
put-1  lne< “power  of  the  'X‘  in  fractional  fore."); 
put_l  ine<" I f  a  sign  is  entered,  it  Must  be  the  "); 
put-i  ine<“ first  character.  Ho  blanks  are  ai lowed.") 
put_l  ineCThe  max  allowable  size  is  Q  digits  per"); 
put-1  ine("nueber ." ); 

put_l  ine(’Sawp!e  entries:  1/2  ,  +1/2  ,  or  -1/2") 
new- 1  ine; 


TO  :■  POUER; 
case  STRUCTURE  is 


when 

T 

»> 

FROTI 

M 

0; 

BY 

m 

i; 

when 

■2' 

■> 

FROM 

m 

o; 

BY 

m 

2; 

when 

'3‘ 

■> 

FROTI 

m 

i; 

BY 

M 

2; 

when  others*> 

FROM 

m 

0; 

BY 

m 

i; 

end  case; 


while  (FROTI  <■  TO)  loop 
G£T_POUER<FROff,  COEFF); 
HRCLRURIH(FROH)  COEFF; 
FROTI  :■  FROTI  ♦  BY; 
end  loop; 


end  GETJCOEFFICIENTS; 


function  PROOUCT(FROM,  TO,  BV:  integer)  return  LONG-FLOflT  is 

RESULT:  LONG-FLORT  I.O; 

LOOP-TEST :  integer  :•  FR0T1; 


begin 

shite  (LOOP-TEST  <■  TO)  loop 

RESULT  :■  RESULT  *  LOMG_FLOflT( LOOP-TEST); 
LOOP-TEST  LOOP-TEST  +  0V; 
end  loop; 

re turn (RESULT); 

end  PRODUCT; 


■  A  .Ik  r. 


-A. 


function  FACTORIAL  (HUMBER:  integer)  return  LONG-FLOAT 
RESULT:  LONG-FLOAT  :«  1.0; 
begin 

for  I  in  2. . NUMBER  loop 

RESULT  :■  RESULT  *  L0NG_FL0AT< I ); 
end  loop; 

re turn< RESULT); 

end  FACTORIAL; 

tnd  CuTiloli-rRoCS; 


—  Date: 

—  Uersion: 

—  Name: 

—  Module  Number: 

—  Description: 


—  Passed  Uariables: 

—  Returns: 

—  Globals  Used: 

—  Modules  Cal  led: 

—  Author: 

—  History: 


28  Noueeber  1985 

1.0 

FUNCTION-PACKAGE 

4.0 

This  package  contains  modules  that  are  called  to  either 
compute  a  predefined  power  series  expansion  of  a  function 
or  allow  a  users  to  enter  their  own 
N/fl 

N/fl  _ 

GLOBAL-DATABASE 

Hone 

Copt.  Steven  R.  Hotchkiss  and 

Copt.  Jennifer  Fried 

0e«Mloped  as  a  thesis  and  ROA  project 


package  FUNCTION-PACKAGE  is 
procedure  SIN-SERIES; 
procedure  TAN-SERIES; 
procedure  RS IN-SERIES, 
procedure  ATAN-SERIES; 
procedure  EXP-SERIES; 
procedure  BUILD-SERIES; 


end  FUNCTION-PACKAGE; 
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GLQBRL-DfiTfiBRSE;  use  GLOBRL-DflTRBflSE; 

COTtlOH-PROCS;  un  COMMON-PROCS; 

TEXT-10;  un  TEXT_I0; 

package  body  FUNCTION-PRCKRGE  is 


procedure  SltUSERIES  is 
N:  integer; 
begin 

—  get  the  posers  of  the  numerator  and  denominator  polynomials. 

—  Also  prompt  the  user  for  a  convergent  epsi Ion. 

POMER-PRONPT<M,K,EPS>; 

—  Compute  the  pomer  of  the  haclaurin  series.  It  is  the  sum  of  the  poser 

—  of  the  numerator,  denominator,  and  the  value  tso 
N  :•  M  +  K  +  2; 

—  Compute  the  initial  approximating  polynomial 
for  I  in  0. .25  loop 

MflCLRURINCI  >  :*  0.0; 
end  loop; 

for  I  in  1 . . (<N+1 )/2>  loop 

flfiCLflUR  INC  1*2-1 )  :«  -1 .0**! /FRCT0RIRL<2*I-1 ); 
end  loop; 


end  SIN-SERIES; 


proc«dur«  TfifLSERIES  is 


N:  integer; 
begin 

—  get  the  posers  of  the  numerator  and  denoeinator  polynomials. 

—  Also  prompt  the  user  for  a  convergent  epsi Ion. 

POUERJPROnPTth,  K,  EPS  >; 

—  Compute  the  poser  of  the  hoc  lour  in  series.  It  is  the  sum  of  the  poser 
— *  of  the  numerator,  denominator,  and  the  value  tso 

N  :■  h  +  K  ♦  2; 

—  Compute  the  initial  approximating  polynomial 
for  I  in  0. .25  loop 

rnCLflURINU)  :>  0.0; 
end  loop; 

MfiCLflURIfKI)  :■  1.0; 

for  I  in  1 . . <<N+1 )/2)  loop 

MRCLRUR I N(2* I + 1 )  :«  PR00UCT<2,  2*1,  2>  / 

FACTOR  I RL<2*I+1>; 

end  loop; 
and  TflhjSERIES; 


procedure  ASINjSERIES  is 


N:  integer; 
begin 

—  get  the  powers  of  the  nueerotor  and  denoeinator  polynomials. 

—  Also  prompt  the  user  for  a  convergent  epsilon. 

POWER_PRONPT<N,K,EPS); 

—  Compute  the  power  of  the  Moclaurin  series.  It  >s  ine  sum  of  the  power 

—  of  the  numerator,  denominator,  and  the  value  two 
N  N  +  K  +  2; 

—  Compute  the  initial  approximating  polynomial 
for  I  in  0. .25  loop 

mCLAURINO)  :»  0.0; 
end  loop; 

for  I  in  1..<<N+1>/2)  loop 

MACLAUR I M< I *2- 1 )  :*  PR00UCT<1,  << l-2)*2+1 >,2)  / 

PA00UCT<2,<l*2-2>,2>  * 

L0HG-FL0ATC 1*2-1 >; 

end  loop; 


end  ASINJSERIES; 


procedure  ATANjSERIES  is 


N:  integer; 
begin 

—  get  the  powers  of  the  numerator  and  denoeinator  polynomials. 

—  Also  prompt  the  user  for  a  convergent  epsilon. 

POUEAJ»ROflPT<M,K,EPS>; 

—  Compute  the  poser  of  the  haclaurin  series.  It  is  the  sue  of  the  poser 

—  of  the  numerator,  denominator,  and  the  value  tso 
M  :■  H  +  K  ♦  2; 

—  Compute  the  initial  approximating  polynomial 
for  I  in  0. .25  loop 

mCLAURIMI)  :*  0.0; 
end  loop; 

for  I  in  1. . <<N+1>/2>  loop 

MACLAURIN<2*I-1 )  :=  -1 .0**< 1-1 > /FACTORIAL <2*1-1 >; 
end  loop; 


end  ATANJSERIES: 


proctdurt  BUI LOOSER I ES  is 

H  :  integer; 

STRUCTURE  :  character; 

begin 

set_page_l  ength<24); 

—  get  the  powers  of  the  numerator  and  denoainator  polynoaials. 

—  Also  proept  the  user  for  a  convergent  epsi Ion. 

POUER-PROMPTOI,  K,  EPS ); 

—  Coepute  the  power  of  the  ttoclaurin  series.  It  Is  the  sue  of  the  power 

—  of  the  numerator,  denoainator,  and  the  value  two 
N  :■  n  +  K  +  2; 

—  Coepute  the  initial  approxiaating  polynoaial 
for  I  in  0. .29  loop 

HRCLRURIfKI)  :■  0.0; 
and  loop; 

—  Proept  the  user  for  the  structure  of  the  polynoaial 
new-page; 

LI: 

loop 

put-1 1 ne< "Enter  1  if  al I  powers  of  X'); 
put-1 i net "Enter  2  if  only  even  powers  of  X">; 
put< "Enter  3  if  only  odd  powers  of  X  — >  "); 
get  (STRUCTURE); 
new-line; 

if  T  >  STRUCTURE  or  STRUCTURE  >'3’  then 
put-1  ine< "Bad  Entry.  Try  again."); 

else 

0ETJC0EFF I C I EHTS< STRUCTURE,  ID; 
end  if; 

exit  LI; 

end  loop  LI; 


end  BUILD-SERIES; 
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procedure  EXPjSERIES  is 


N:  integer; 
begin 

—  get  the  posers  of  the  numerator  and  denoeinator  polynomials. 

—  Also  proept  the  user  for  a  convergent  epsilon. 

POHER_PROHPT<H,K,EP$); 

—  Coepute  the  poser  of  the  hoc  I  our in  series.  It  is  the  sue  of  the  poser 

—  of  the  numerator,  denoeinator,  and  the  value  teo 
N  :■  H  +  K  +  2; 

—  Coepute  the  initial  approxieating  polynoeial 
for  I  in  0. .25  loop 
riRCLfURlhO)  0.0; 
end  loop; 

for  I  in  0. .N  loop 

nflCLflURINCI)  :»  1.0  t  FftCTOR I «.<!>; 

end  loop; 

end  EXPjSERIES; 
end  FUNCT 1 ON-PRCKRGE : 


—Date:  28  Noveeber  1985 

—  Version:  1.0 

—  Naee:  APPROXIMATORS 

—  Module  Muaber:  5.0 

—  Description:  This  package  i  ncludes  procedures  that  coepute 

—  approx i eat  ions  to  user  selected  functions 

—  Passed  Variables:  N/fl 

—  Returns:  M/fl 

—  Global*  Used:  GL08AI _ DfiTRBRSE 

—  Modules  Cal  led:  N/fl 

—  Author:  Capt.  Steven  A.  Hotchkiss  and 

—  Capt  Jennifer  Fried 

—  History:  Developed  as  a  thesis  and  AOfl  project 


package  APPROXIMATORS  is 

procedure  COMPUTE-TCHEBYSHEU ; 
procedure  MENUtCHOICE:  out  character); 
procedure  C0MPUTE_PA0E_RPPR0XIHATI0NS; 
procedure  COMPUTEXMK; 


and  APPROXIMATORS; 


•ith  GLOBAL-DATABASE;  um  GLOBAL-DATABASE; 

•ith  FUMCTIOH-PACKAGE;  um  FUNCTION-PACKAGE; 

•  ith  TEXT-10;  um  TEXT-10; 

package  body  APPROXIMATORS  is 

package  FLT-IO  is  now  FLOAT- 1 0 < LONG-FLOAT ) ;  um  FLT-IO; 

procedure  CONPUTE-TCHEBYSHEU  Is 
begin 

—  build  the  global  table  "T“  containing  the  coefficients  for 

—  each  of  a  series  of  Tchebyshev  polynomials 

T<0,0>  :■  1.0; 

T< 1, 1 >  :•  1.0; 

T<2,0>  ;■  -1.0; 

T<2,2)  :■  2.0; 

for  I  in  3. .25  loop 

for  J  in  0. .25  loop 

T<1 , J>  :«  T<I,J>  -  T<I-2,J>; 
end  loop; 

for  J  in  0.  .24  loop 

T< I ,  J+1  >  :■  T< I ,  J+1 )  +  2.0  *  T< 1-1, J>; 
end  loop; 

end  loop; 

end  COHPUTE-TCHEBVSHEU; 
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procedure  MEM) (CHOICE:  out  character )  is 

OUT-CHOICE:  character; 

BflOJCHOICE:  exception; 

begin 

set_page_l  ength<24  >; 

—  clear  screen  and  print  eenu 
nee-page; 

put-1  ine< “Choose  function  to  be  approximated” >; 
nee_l  ine; 

put_l  ine< “Enter  1  for  sin”); 

put-1  ine( “Enter  2  for  tan’); 

put_  I  ine  ("Enter  3  for  arcsin“>; 

put-1  ine< “Enter  4  for  arctan“>; 

put-1  ine(  “Enter  5  for  exp“>; 

put_l  ine< “Enter  6  for  user  defined  function” > 

put-1 ine( "Enter  7  to  quit*); 

loop 

for  I  in  MRCLRUR IN’ range  loop 
MRCLflURIMI)  :■  0.0; 
end  loop; 

nee_l ine; 
put(“— >  “); 
get(OUTjCHOICE); 

CHOICE  :»  OUT-CHOICE; 
case  OUT -CHOICE  is 

ehen  ' I*  ■>  SIN-SERIES; 

•hen  ‘2‘  ->  TRN-SERIE8; 

•hen  *3'  »  RSI  MISERIES; 

•hen  ‘4'  ->  ATflNjSERIES; 

•hen  ‘3‘  ■>  EXP -SERIES; 
when  ‘O'  ■>  BUILD-SERIES; 

■hen  ■>  null; 

■hen  others  ■>  raise  BflOJCHOICE; 

end  case; 

exit; 

end  loop; 
exception 

■hen  BflOJCHOICE  » 

put-1 ine(“lnual id  entru.  Try  again*); 


procedure  COrfH/TE_PfCE_flPPROXIMRTIOre  is 


NJ1RX:  integer; 

NUM  :  integer  :>  0; 

OEM  :  integer  :»  t; 

TEMP  :  L0NG_FL0flT; 

UORK  :  MATRIX; 

B  :  UECTOR; 

begin 

—  this  procedure  converts  the  initial  approximating  polynomial 

—  (the  Mac  I  our  in  pay  series)  into  a  rational  approximation 

—  clear  out  this  PROS  approximation's  numerator 

—  and  denominator  polynomials 
for  SERIES  in  0. ,25  loop 

for  NUIUEN  in  0..1  loop 

for  COEFFICIENT  in  0..25  loop 

R(S£RIES,NUN_DEN, COEFFICIENT)  :■  0.0; 
end  loop; 
end  loop; 
end  loop; 


for  I  in  0..M  loop  —  loop  for  all  posers  of  the  numerator 

for  J  in  0.  K  loop  —  loop  for  all  powers  of  the  denominator 
If  (I  >■  J)  then 

—  build  a  work  matrix  to  solve  simultaneous  equations 
B<0)  :«  1.0; 

NJIflX  :«  I  ♦  J; 

for  S  in  0. . (N-MAX  -  I  -  1)  loop 
for  N1  in  0. .J  loop 

U0RK(S+1,N1 )  :»  MRCLAUR I N(abs(N_MRX  -  S  -  N1)); 
if  <N1  -  0)  then 

B(Stl)  -MflCLRURIN<abs<NJ1flX  -  S  -  Nl)>; 
end  if; 
end  loop; 
end  loop; 

—  Solve  simultaneous  equations  for  denominator  coefficients 
for  Ml  in  1 . . J  loop 

if  <W0RK(N1,N1>  -  0.0)  then 
SETUP: 

for  N2  in  1.  .J  loop 

if  <U0RK(M2,N1)/-  0.0)  then 
TEMP  :•  B(N2); 

B(N2)  B(N1); 

B(N1)  :■  TEMP; 
for  N3  in  1. .  J  loop 
TEMP  :•  M0RK(M2,M3>; 

WORK < M2, M3)  :•  U0RK(N1,N3); 

W0RK(N1,N3>  TEMP; 
end  loop; 
exit  SETUP; 
end  if; 

end  loop  SETUP; 


•nd  if; 


i  %' 


TEMP  :«  U0RK(N1,  N1); 
if  TEMP  /■  0.0  then 
B<N1)  :■  B(M1)/TEMP; 

•  Is* 

B<N1>  :■  0.0; 

•nd  If; 

for  M2  in  1. .J  loop 
if  TEMP  /•  0.0  than 

W0RK<M1,N2>  :■  U0RK<M1,N2)/TEHP; 

•Is* 

U0RK(N1,N2)  :■  0.0; 

•nd  if; 

•nd  loop; 

for  N2  in  1. .J  loop 
if  <N1  /»  N2>  then 
TEMP  :»  -U0RK(N2,N1 >; 
for  M3  in  1. .J  loop 

U0RK(N2,N3)  :»  U0RK(N2,M3>  +  U0RK<N1,N3>  *  TEMP; 

B<M2)°^'b<M2)  +  B<N1>  *  TEMP; 

•nd  if; 

•nd  loop; 

•nd  loop; 

—  us*  denominator  coefficients  to  compute  the  numerator 

—  coofficiants,  and  build  th«  series  of  PRDE  approx i eat  ions 
for  Ml  in  0. . I  loop 

for  N2  in  0. .N1  loop 

R<I+J,MUM,M1>  R< l+J,MUM,M1 >  +  B<N2>  *  MRCLRURIN(M1-M2>/ 
B(0>; 

•nd  loop; 

•nd  loop; 

for  Ml  in  r*v«rs«  0..J  loop 
R(I+J,0EM,N1>  :•  B(Nt)/B<0>; 

B(M1 )  :■  B(M1 )  /  6(0); 
end  loop; 

—  Coeput*  th«  D’s  that  am  used  to  coeput*  C(»,k) 

—  0<I+J+1)  ■  SUMIL=0  to  J  (Mac  I  our i n< I +J+ 1-L  >*B<L  > 1 
D< I +J+ 1 )  :*  0.0; 

for  L  in  0. .J  loop 

0<I+J+1>  D(l+J+1>  +  MRCLRURIMC I+J+1-L)  *  B(L>; 

end  loop; 


•nd  if; 

•nd  loop; 

•nd  loop; 

end  COMPUTE-PROE-flPPROXIMRTIOMS; 


/'•  /■  .*•  -  ■  .  ■  ,>  .**  .*■ 
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proctdurt  COMPUTE-CMK  is 


A:  integer  :■  0; 
B:  integer  :»  1; 
LRMOR:  UECTOR; 

begin 


—  Cowxtta  the  Laedas  (a  I  oho*  1) 

LRHOR<0>  :■  -<0<M+K+1>  *  T<M+K,0>>  /  <2.0**<M+K)>; 

for  J  in  0. .<M+K-1 >  loop 
if  D<J+1>  /■  0.0  then 

LRM0fl<J+1>  <0<H+K+1>  *  T<M+K+t,J+1>>  t  <<2.0  **<M+K>>  *  0<J+1>> 

else 

LRMDR<J+1>  :«  0.0; 
end  if; 
end  loop; 


—  Load  Pe<X>  and  Qe<X)  eith  their  R  and  B  coefficients  respectively 
for  I  in  0.  .11  loop 

C<fl, I)  :•  R<M+K,fl,l>; 
end  loop; 

for  I  in  0. .K  loop 

C<B, I )  :■  R<H+K,BJI); 
end  loop; 

—  Coepute  coefficients  “fl*  of  nueerator  and  “B"  of  denoeinator 
for  J  in  0..O1+K-1)  loop 

for  K  in  0. .25  looo 

R<J,R,K>  R<J,fl,K>  *  LRM0fl<J+1>; 

C<fl,K>  :■  C<fl,K)  ♦  R<J,fl,K>; 

R<J,B,K>  :■  R<J,B,K)  *  LRM0fl<J+1>; 

C(B,K)  :*  C(B,K>  ♦  R<J,B,K>; 
end  loop; 
end  loop; 

C<R,0)  :■  C<fl,0>  +  LRHOR(O); 

for  I  in  reverse  0. .25  loop; 

C<fl,  I)  :»  C<fl,  IVC<B,t»; 

C<B, I )  :■  C<B,l)/C<B,0>; 

end  loop; 


end  COMPUTEJOK; 
end  APPROXIMATORS; 
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